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Abstract 

We consider the bilocal conductivity tensor, the two-probe conductance and 
its fluctuations for a disordered phase-coherent two-dimensional system of 
non- interacting electrons in the presence of a magnetic field, including cor- 
rectly the edge effects. Analytical results are obtained by perturbation theory 
in the limit a xx ^> 1. For mesoscopic systems the conduction process is dom- 
inated by diffusion but we show that, due to the lack of time-reversal sym- 
metry, the boundary condition for diffusion is altered at the reflecting edges. 
Instead of the usual condition, that the derivative along the direction normal 
to the wall of the diffusing variable vanishes, the derivative at the Hall angle 
to the normal vanishes. We demonstrate the origin of this boundary condition 
from different starting points, using (i) a simplified Chalker-Coddington net- 
work model, (ii) the standard diagrammatic perturbation expansion, and (iii) 
the nonlinear <r-model with the topological term, thus establishing connec- 
tions between the different approaches. Further boundary effects are found 
in quantum interference phenomena. We evaluate the mean bilocal conduc- 
tivity tensor a^ u (r,r'), and the mean and variance of the conductance, to 
leading order in l/o~ xx and to order (cr xy /a xx ) 2 , and find that the variance 
of the conductance increases with the Hall ratio. Thus the conductance fluc- 
tuations are no longer simply described by the unitary universality class of 
the o~ X y — case, but instead there is a one-parameter family of probability 
distributions. Our results differ from previous calculations, which neglected 
(Tajy-dependent vertices beyond the change in boundary conditions. In the 
quasi-one-dimensional limit, the usual universal result for the conductance 
fluctuations of the unitary ensemble is recovered, in contrast to results of pre- 
vious authors. We also give an extensive discussion of current conservation 



* Present Address: International Center for Theoretical Physics, 34100 Trieste, Italy 



1 



conditions in the nonlinear sigma model. 
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I. INTRODUCTION 



In the past decade much attention has been given to the statistical properties of quantum 
conductors with complete phase coherence (with size L smaller than the phase coherence 
length Li n ). A notable feature of such systems is the lack of self-averaging in their transport 
properties. In mesoscopic systems for which L is less than the localization length £, the 
conductance fluctuation amplitude (the standard deviation of the conductance), is found, 
at low magnetic field, to be of order 1 and independent of system size and the degree of 
disorder, (but dependent on the dimensionality, shape, and overall symmetry properties of 
the system) |TJ-0. (Note that, in the present paper, all conductivities a and conductances g 
have a factor e 2 /h removed, so they are dimensionless in two dimensions; their dimensionful 
analogues are recovered by multiplying by e 2 /h.) These universal conductance fluctuations 
(UCF) have been well understood within the framework of perturbation theory ]3],|],[7] and 
the one parameter scaling theory of quantum conductance The physics underlying 

the UCF is the long-ranged spatial correlation among the wave functions of the conduction 
electrons in the diffusive regime. The universality of the phenomenon has also stimulated 
the formulation of a random matrix theory description of quasi-one- dimensional (quasi-lD) 
|T0Hl3| and quasi- zero-dimensional (quasi-OD) conductors []14],jl5], which reproduces quanti- 
tatively the results of diagrammatic perturbation theory. In the earlier perturbative work 
the effect of a magnetic field was only included through the introduction of an appropriate 
Aharonov-Bohm phase in the zero-field propagator. This leads to the elimination of the 
Cooperon contributions to the conductance fluctuations, causing a cross-over from the so- 
called orthogonal to the unitary ensemble, and a consequent reduction of the variance by a 
factor of two |||7|]. In the quasi-lD case, this result is also recovered by the random-matrix 
approaches |T^JT3|] . 

In two dimensions (2D), the interplay of quantum interference effects and magnetic field 



leads to the quantum Hall effect in high magnetic field B [fLq| . In mesoscopic samples, 
conductance fluctuations persist in fairly high magnetic field u c tq > 1, where u c is the 
cyclotron frequency, tq is the elastic scattering time in zero magnetic field, and the fluctuation 
amplitude remains comparable to the low field limit ]T7|-P^]. For B fields sufficiently high 
that quantization of the Hall conductance sets in, the fluctuations are strongly suppressed in 
the plateau regions, but reappear in both longitudinal and Hall resistance in the transition 
regions between plateaus. It is therefore of theoretical interest to generalize the theory 
for conductance fluctuations to all fields. At lo c tq > 1, the trajectories of the electrons 
are significantly influenced by the Lorentz force between successive scattering events and 
the dynamical effect of the magnetic field must be treated. The diffusion at high field 
occurs by a different mechanism from the low field regime. For the short-ranged random 
potential model, the center of the cyclotron orbit hops a length of order the cyclotron radius 
R c whenever it encounters a scattering center, therefore R c plays the role of the mean 
free path I. The bare conductivity a® x in the middle of the iVtli Landau level is of order 
(2N + l)/n (N — 0, 1, . . . ). Despite the altered nature of the microscopic diffusion, a 
unified treatment of mesoscopic conductance fluctuations in relatively high fields is possible 
because, even at uj c tq > 1, there exists a perturbative regime where the transport process is 
dominated by diffusion. As long as there are many Landau levels occupied, l/cr®x or 
serves as a small parameter, and perturbation theory is still useful. Previous perturbation 
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theories have shown that the weak localization correction to a xx for the unitary class is of 
order — (a® x )~ 1 log(L/l) P5-271- The localization correction is relatively small for systems 
with L less than the crossover length £ pert = le^ xx ' (for a® x large). For I < L < £ pert , 
the conductance fluctuations are expected to be similar to the UCF. At L > £ peri , the 
renormalization group flows for the system, driven by non-perturbative effects [p8|-|3T]j, carry 
it either to one of the localization fixed points where a xx — > and a xy becomes quantized, or 
to one of the non-trivial fixed points where a xy is half integer and a xx approaches a universal 
value. Numerical work finds that the conductance fluctuations in the critical regimes [p0| -[23 



have a different distribution from the UCF and one would expect this distribution to be 
beyond the scope of a perturbative treatment. We comment further on the critical regime 
in the conclusion, Sec. VI. 

In this paper we study the conductance and its variance for a two-probe geometry in two 
dimensions in the presence of disorder and magnetic field and in the perturbative regime 
(/ < L < ipert)- As shown in Figure |I|, the edges of the sample are defined by a hard- wall 
confinement potential and the two ends of the sample are connected to highly conducting 
leads. The first analytic work on this problem was by two of the authors (Xiong and Stone 
1 32[) in which they generalized the previous diagrammatic perturbative techniques |27j to 
treat the conductance fluctuations. At the level of the self-consistent Born approximation 
(SCBA), the only effect they found of the magnetic field B was a field-dependent diffu- 
sion constant D(B). Since the value of the diffusion constant cancels from the conductance 
fluctuations, they found no effect of the magnetic field on the amplitude of the two-probe 
conductance fluctuations (other than the well-known factor of two reduction associated with 
the crossover to the unitary ensemble), although the correlation field, B c , which determines 
the spacing of the fluctuations in magnetic field, was found to increase with increasing field 
in a manner consistent with experiment [[Tj|. The reason for the increase is that for systems 
with L > L in , B c ~ (fi /Ll n , where L 2 in = D(B)ri n (<po = h/e, r in is the inelastic scattering 
time). Since the diffusion constant D(B) decreases with increasing magnetic field (a result 
reviewed in Sec. II), B c increases. Although this conclusion about the correlation field is ba- 
sically correct, the conclusion concerning the fluctuation amplitude is now understood to be 
correct only for a periodic boundary condition in the transverse direction and must be revised 
for the case of a system with reflecting edges. As discovered independently by Khmel'nitskii 
and Yosefin (KY), Maslov and Loss (ML), and one of the present authors (Read) f33|-|35 



the boundary condition is modified from the vanishing of the normal derivative of the dif- 
fusing variable to the vanishing of the derivative at an angle to the normal. KY and Read 
further showed that this angle is the Hall angle 9jj = tan -1 cr° y /cr xx , where a xy is the bare 
Hall conductivity. These authors pointed out the possibility that the mesoscopic conduc- 
tance fluctuations may depend on magnetic field due to the boundary condition. KY and 
ML attempted to evaluate this dependence both numerically and analytically. Simulations 



performed for the two-probe conductance of small systems in the non-quantized regime \M 
show that the maximum fluctuation amplitude appears towards the bottom of the Landau 
levels where the Hall ratio is large, indicating some dependence on the Hall ratio. How- 
ever, the analytic calculations by KY and ML p3|^4| do not agree with our present results 
since these authors have merely modified the diffusion propagator in previous expressions 
for UCF diagrams. Like KY, we find that a xy enters not only the boundary condition but 
also the current vertex. Moreover the altered boundary condition permits new diagrams to 
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occur which roughly speaking describe interference effects associated with a a xy dependent 



"interaction" of the diffusion modes. These diagrams, which were not considered by ML 
and KY, must be included when edges are present. We evaluate the two-probe conductance 
and its variance to leading order in l/c xx and to order 7 2 , where 7 = cr xy /a xx . For wide 
samples with W ~ L, where L and W are now the length and the width of the sample, we 
find that the variance does not depend on a xx and a xy individually, but does depend on 7. 
The variance increases as 7 2 for small 7 and hence is no longer independent of magnetic field 
(although it is still independent of size in this order in l/a® x , and has no direct dependence 
on the mean- free path). Interestingly, however, in the quasi-lD limit (W <C L), the Hall 
ratio is absorbed into an effective ID conductivity which cancels in all diagrams. Therefore 
in this limit the UCF result of the unitary class is recovered (in contrast to the claims by KY 
and ML that this is modified for 7 7^ 0). The implication is that quasi-lD conductance fluc- 
tuations are still described by the standard random matrix theory of disordered conductors, 
but that the 2D fluctuations, even in perturbation theory, define a family of random-matrix 
ensembles parametrized by the Hall ratio 7. 

In this article, we use the disorder-averaged diagrammatic approach and the field theo- 
retical approach in a complementary way. To study transport properties of a system with 
phase coherence, the appropriate starting point is the bilocal conductivity tensor a^r, r'). 
In Sec. II, we evaluate the mean o" MV (r, r') to leading order in l/u xx using the diagrammatic 
approach and demonstrate the microscopic origin of the edge contributions. In Sec. Ill, we 
set up the field-theoretical formalism for the evaluation of linear response functions. We 
discuss the connection between the tilted boundary condition and the nonlinear cr-model 
action with a topological term proportional to a xy p8| , |30[| . Previously it was known that 



this topological term is crucial to the critical transition of the quantum Hall effect at large 
length scale (L ^> ^ pert ) [ |2T||50| ] . For a system with reflecting edges, this term is a non- 



vanishing surface term, which does influence transport properties in perturbation theory, 
valid when L <C (, per t, through various cr^-dependent boundary contributions. In Sec. IV 
and Sec. V, the conductance and its variance are calculated by expanding in power series 
in 1/(T XX and 7. In the remainder of the introduction we begin the discussion of the main 
ideas and summarize our results. The details of the calculations, and further discussions, 
are given in later sections. 



A. Local conductivity parameters, two-probe conductance, and edge states 

In this article we focus on calculations of the two-probe conductance in a magnetic field. 
It is however possible to generalize our calculations to treat the conductance matrix of a 
multi-probe conductor as has been done previously for zero (or weak) field p6||37| . Two- 



probe conductance describes an experimental set-up in which voltage measurements are 
made only between the current source and sink and not between distinct voltage probes as 
in a typical Hall measurement (two-probe measurements are not uncommon for mesoscopic 
conductors, because of the difficulty of making multiple contacts). Therefore such a mea- 
surement cannot separately determine a xx {B) and a xy (B). In this subsection, we show how 
the assumption of a local form for the conductivity in a system with edges leads to the result 
that the two probe conductance is approximately proportional to a xx when a xx ^> \cr xy \ and 
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to \(T X y\ when \a xy \ 3> a xx . Interestingly, such an argument already indicates the appearance 
of the tilted boundary condition which affects the conductance fluctuations as well. 

We wish to find the current produced in linear response to an applied electric field. We 
will consider the two-probe conductance which results from assuming a local form for the 
conductivity: 



e 2 



j» = -a» v E v (r), (1.1) 

where due to the macroscopic homogeneity of the sample and Onsager relations the con- 
ductivity parameters obey a^x = a yy an d °~x y = ~ a yx- ft * s a common misconception that 
£ in this formula is the same as the applied electric field E. In fact, in general, £ in this 
formula should be interpeted as the electromotive field, that is, as minus the gradient of the 
electrochemical potential V = — /i/e, where is the electric potential, and /x the chemical 
potential (more generally, S = E + V/x/e). In this paper, we will neglect electron-electron 
interactions of all kinds, so <fi can be viewed as the externally-applied electric potential. (If 
it is desired to include Coulomb interactions through a self-consistent potential, then <fi is 
the total electric potential, the sum of the externally-applied potential and the potential 
produced by the electrons in the system; the latter potential is determined by the change in 
the expectation value of the density, in response to the external field, through the 3D Pois- 
son equation. This is distinct from similar-looking equations below which have the form of 
the 2D Laplace equation. Other short-range interaction effects would contribute to fi.) The 
chemical potential in the above expressions is defined in terms of a local quasi-equilibrium 
(in the presence of a nonzero response current) which must be established by inelastic effects. 
Hence, this formulation is only valid on scales greater than the inelastic length L in (this of 
course requires that there be some interaction between the electrons). It is only in this sense, 
which implies the absence of effects due to single-particle phase coherence, that Eq. (|1.1|) is 
a classical formula; it does not require that all effects of quantum mechanics be neglected. 
Eq. (|1 . 1|) (when valid) is the most convenient form for expressing the linear response, since 
a voltage measurement determines electrochemical potential differences. It does not imply 
that a local relation exists between the electromotive and electric fields. As the chemical 
potential is determined by the local conditions, in particular by the local density, and that 
density is affected by the transport, which in our case will be diffusive, this relation is not 



local. Thus the current response to the electric field is actually nonlocal, as in Ref. p8 
even in this "classical" case. We will return to this in Sec. ID. 

We will now calculate the two-probe conductance of a rectangular sample with insulating 
edges connected to conducting leads at each end. The potential in the leads is assumed to 
be held at constant values, V\ and Vi. At every point on the insulating edges the normal 
current must vanish. Using Eq. ( |L1|) and S = — VV it follows that (d y — jd x )V(r) = 0, i.e. 
the electromotive field at the Hall angle to the normal must vanish. Thus in this case the 
appearance of a tilted boundary condition on the electrochemical potential follows simply 
from the fact that the field is tilted from the current by the Hall angle. From the continuity 
equation, V • j = 0, and one finds that V 2 V = in bulk. Solving the Laplace equation for 
V(r) with fixed voltages at the two ends and the tilted boundary conditions at the edges 
is not a simple exercise, but has been done for this 2D rectangular geometry by conformal 
mapping PUJ (we give a solution in another form in Sec. ID). From this solution one can 
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obtain the two-probe conductance for arbitrary Hall ratio, however here we analyze only 
its limiting behavior. The conductance g° can be found by integrating the current over a 
cross-section perpendicular to the current flow: 



g° = -a° xx / dy (d x + 7^) V(r)/(V 2 - V,). (1.2) 

First consider the case 7 — > 0; it is then convenient to integrate over all transverse cross- 
sections in the sample and divide by the sample length L. The integral then just evaluates 
the voltage at the two ends where it is fixed, yielding the familiar ohmic result 

\hag = o°„W/L. (1.3) 

7^0 



In the opposite limit a xx — > 0, \j\ — > 00, the boundary condition implies that everywhere 
along the edge S x = 0, i.e. there is no voltage drop along the edge and the voltage must 
remain equal to V\ along one edge and V 2 at the other edge except at singularities at 
diagonally opposite corners (where it jumps between V\ and V 2 ). The transverse potential 
drop is equal to (Vi — V2)sgn7, therefore 



rW 

lim g° = - lim a° xx / d yi d y V(r)/(V 2 - V,) = \a° xy \. (1.4) 



So, using the classical local conductivity, the two-probe conductance changes from being 
dominated by the longitudinal conductance at small 7 to being dominated by the Hall 
conductance at large 7. We can also solve the limit L/W — > 00 with 7 fixed; for a fixed 
current, the voltage drop is dominated by the part of the geometry a distance greater than 
W from the ends, in which the current distribution is essentially independent of y, and one 

finds that g° -> [{(<r°J 2 + ( a l y ) 2 } I a lx] w I L = ( l /plx) w / L - Note that the crossover to this 
behavior occurs at a value of L/W that depends on 7. In contrast to the above results for a 
rectangular sample with edges, for periodic boundary conditions in the transverse direction, 
which is equivalent to transport along a cylinder, the conductance is always a xx W/L, for 
any value of 7. We note that this geometry is equivalent, through a conformal mapping, 
to the Corbino disk geometry, in which the voltage drop is radial, and since the equations 
are conformally invariant, results for the cylinder apply to the disk also. Although the local 
formulation cannot be used in the fully phase-coherent ("quantum") case where L in ^> L, we 
expect that the physics illustrated by this argument will be relevant to the average quantum 
conductance. 

Indeed, in the quantized Hall regime L>(, previous arguments based on the Landauer 
formula for two-probe conductance in terms of transmission coefficients have noted the 
relation between two-probe conductance and Hall conductance. These approaches assume 
that the incident and outgoing channels are iV edge states [^0| , p| . These edge states are 



analogous to classical skipping orbits advancing in one sense along each edge and occasionally 
being scattered into the bulk. Any actual calculation of these transmission coefficients will be 
equivalent to evaluating the bilocal conductivity between the two ends [i2|, however physical 



arguments are made that in high field the backscattering of edge states will be suppressed 
giving perfect edge transmission and g = N = \cr xy \; i.e. the two-probe conductance is equal 
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to the Hall conductance which takes its quantized value. (Since a xx = this is consistent 
with the classical result above). Although physically appealing, such arguments assume 
that the interaction of bulk and edge states can be ignored. However in a bulk 2D sample 
it is known that it is only localization effects which prevent edge states from backscattering 
through the bulk states; they thus require that W 3> £, the localization length. We note 
that away from the critical values E c n (N = 0, 1, . . . ) of the Fermi energy, which lie near 
the center of the iV-th Landau level, £ ~ £, per t, while £ — > oo as Ep — > E c ^ for each N. We 
cannot calculate analytically the conductance in this non-perturbative regime L, W > £, per t', 
however our results below do describe samples for which L, W <C £ pert and localization 
effects have not inhibited backscattering of edge states. The key ingredient to describe the 
edge-bulk coupling in the perturbative regime is the tilted boundary condition for diffusion 
which we derive in the next subsections. In the fully phase-coherent case, the boundary 
condition is however modified further. 



B. Classical network model, edge states, and tilted boundary condition on diffusion 



Khmel'nitskii and Yosefin J53J (KY), Maslov and Loss |34] (ML), and Read [35[ have 



previously obtained the boundary condition on the diffusion process in high magnetic field. 
KY's argument appears to be related in part to the classical conductivity formulas reviewed 
in Sec. IA. ML considered the effect of the edge in high field on the microscopic diffusion 
process using a Boltzmann equation approach. They found that the tendency to skip in only 
one direction when colliding with the edge does lead to the tilted boundary condition on the 
diffusion equation. They expressed the tilt angle in terms of the ratio of the mean free path 



along the edge and the bulk mean free path. Read [^] used the non-linear sigma model 
approach, which will be described later in this paper. KY and Read were able to identify 
the tilt angle as the Hall angle. In this subsection we rederive the boundary condition in 
a particularly transparent manner using a classical version of the network model for high 
field transport introduced by Chalker and Coddington 13]. In this case one can also see 
immediately that the tilt angle is the Hall angle. 



The original Chalker-Coddington model P3| , fH|| describes the quantum tunneling between 
the semi-classical orbits along the equipotential contours of the smooth random potential 
(see figure 0). To derive the diffusive behavior of the probability density in this model we 
will neglect interference effects and describe each node by the probability that a walker 
approaching it makes a step to the right (R) or left (T); T + R — 1. (This simplified model 
has been considered by several earlier authors |^5)-0; it is essentially classical and could 
serve as a lattice realization of the classical behavior discussed in Sec. IA.) The links of the 
lattice can be divided into four sublattices a = A, B, C, D (see Fig. 2), and each unit 
cell of the lattice contains one of each of the four classes of links. The nearest neighbor 
separation is a. We use p a (i,j,t) to denote the probability of being at link a, site i,j at 
time t. Assuming that it takes time r for a particle to move from one link to the next, one 
can define a random walk problem on the network and write down a probability evolution 
equation: 



PB(i, j, t + r) = p D (i, j, t)T + p A (i, j, t)R, 
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Pc(i,j,t + t) = p A (i,j,t)T + p D (i,j,t)R, 
p A (i,j,t + r) = p B {i,j - l,t)T + p c (i + l,j,t)R, 
pD(i,j,t + r) = p c (i,j + l,t)T + p B (i - l,j,t)R. 



;i.5) 



This problem differs from the usual random walk in the respect that on each link the walk 
is in only one direction, hence breaking time-reversal symmetry. 

The above equation can be diagonalized in the Fourier space. One can show that the 
long-time large-distance modes have a diffusive spectrum 

— itUk = Dk 2 

for uj <C 1/t and k <C 1/a, where D, the diffusion constant, is given by D = ^a 2 RT/ (R 2 + 
T 2 )t. The associated eigenmodes are of the following form in Fourier space: 



(PA(k)\ 


( 


pB(k) 




Pc(k) 




\p D {k)) 





a° xx ik x a + a° xy ik y a 



\ 



1 + 



xy 

-a 



XX 





XI/ 



l)ik x a + o xx ik y a 



cr°y)ik y a 



J 



where in anticipation of our discussion below we identify the two constants a xx = RT / (R 2 + 
T 2 ) and a x = —T 2 /(T 2 + R 2 ) as the the bare longitudinal and Hall conductivity of this 



model. Note that these conductivities satisfy the "semicircle relation" , 

a 2 xx + (a xy + N + 1/2) 2 = 1/4, 



;i.6) 



with N = here in the case of the lowest Landau level, which has been claimed to be a 
general, exact result in the quantum Hall effect [J47|J . In real space, all four components of 
the probability distribution satisfy the same coarse-grained equation: 



- DV 2 p(r,t) = -d t p(r,t). 



;i-7) 



At the absorbing ends, the particle moves away from the tunneling region with constant 
velocity. The fact that the re-entry probability is zero gives the boundary condition at the 
leads: 



p = in the lead. (1.8) 

Since the particles always move in the direction of the arrow on a link, the density on a link 
can also be considered as its current. The differences among the four components define the 
diffusion current densities, which are suitable for coarse graining. For instance, we can define 
j y = [ps(«, j) - PD(i + 1, j)]/a and j x = [p D {h j) - pc{hj)]/ a - For the low-frequency and 
long- wavelength modes, we can show that j y = —(<J xx d y — a xy d x )p(r, t). Along the reflecting 
walls, at zero frequency, the normal current is zero, e.g. psihj) — Pr>(* + l>j) = at top 
edge, which gives the boundary condition: 

(d n - 1 d t )p = 0, (1.9) 
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where n is the outward normal, and t = n x z is the tangential direction of the boundary, 
and 7 = —T/R = (J xy /a xx with our above identification of the bare conductivities in the 
model. These expressions for the current density are of the form used by KY and in Sec. 
IA. 

It is interesting to note that, in the network model, the definition for the local diffusion 
current density is not unique. One can also define, e.g., j' x = [poii, j) — Pc(hj + j' y — 

[psih j) —PD^h j)]/ a - j'y an d j y differ by a total derivative term d x p, and the corresponding x 
components differ by — d y p and by a 5-function boundary term p(x, W)5(y— W) — p(x, 0)5(y). 
In the presence of such an edge current, the boundary condition that ensures "current 
conservation", which in the interior is the requirement V ■ j' = 0, becomes 

r w+o+ 

d x dy'f x (x,y')-f y (x,W) = (1.10) 

Jw-o+ y 



at the top edge, which is still equivalent to ( |1.9| ). These forms for the current density are 



similar to those of ML. The edge contribution ensures that the total current through a 
cross-section transverse to the x-direction is the same, whichever definition of current is 
used. 

One can see that in the long-time, large-distance limit the lack of time-reversal symmetry 
affects the diffusion process only through the tilted boundary condition, which is present 
only because of the edges. As we will show explicitly later, these boundary conditions, 
although derived from a lattice here for convenience, are quite general for conduction with 
broken time-reversal symmetry. As one might hope, in the long-time, large-distance limit 
the microscopic details of different models cease to matter. 

As noted just above, the boundary condition derived in Eq. (1.7) applies for 7 = —T/R, 
the single-node transmission and reflection coefficients of the classical network model. We 
now must further justify identifying this ratio as the bare Hall ratio. There are two ap- 
proaches to this. There is no applied electric field or electric potential in this problem so 
far, so one may simply define the chemical potential on each link as proportional to the 
current (or density) there (in analogy to the Landauer approach for the entire sample). This 
was done by Kucera and Streda and leads to the formulas for a xx and a xy given above. 
In our view a somewhat more satisfactory method is to calculate the steady-state current 
for the network under periodic transverse boundary conditions (i.e., a cylindrical system), 
when current is injected at only one end of the network, with unit current on each incoming 
link at that end. This is just the appropriate time-independent solution of Eq. ( |1.5| ); in the 
absence of edges the solution is the linear k = mode: 

p B (x,y) = hx + bo, 

where b x = l/{L/a + a° xx ), b = 1 - [a° xx /(L/a + a° xx )], and 

p A (x,y) - p B (x,y)\ I ol x b x 

pc(k) - p B (x, y) = (a° xx - a° xy - l)h 
Pd(x,v) - Pb{x,v) ) V (-0-^-1)61 

where the constants cr xx ,a xy are as defined above Eq. (1.6). The above solution has a 
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uniform current distribution with a total longitudinal current, for L 3> a, 

w °L w o 

r = o = — ~ — rr 

and a total transverse current, circulating around the cylinder, 

I a 

*y — u xyi 

justifying the interpretation we have assumed. 

To solve for the current in the presence of edges is much more difficult, even in this 
classical model. We will see, however, that on large scales the problem is equivalent to that 



solved by Rendell and Girvin |}9j and discussed in Sec. IA. The reason for the similarity is 
already clear; in the classical network model, we assumed that the chemical potential was 
proportional to the density, and, when the external electric field is zero, the coarse-grained 
current density was related to this potential in the same way in both cases. 



C. Continuum action and propagator 

The diffusive behavior generated by the classical network model of the previous subsec- 
tion can be reproduced in a simple continuum field theory. Consider the action 

a r ex r 



So = —f- J d'rd^z-^f J d z rt, v d a zd u z. (1.11) 

Here z(x) is a complex scalar field, z is its complex conjugate, and the geometry is the 
same as in Sec. IB. The second term is clearly a total derivative. To obtain the equations 
of motion and boundary conditions in this model, we first rewrite So as 



J d 2 r [zV 2 z] -?fjd 2 r [dMd^ + je^d v )z}}. (1.12) 



Then one can see that the second term is again a total derivative and it can be written as a 
boundary term. Taking the functional derivative with respect to z, we obtain the equation 
of motion for z 

a 

- -^V 2 z(r) = in the bulk 

4 v ' 

(d n — l9t)z(r) = at the reflecting walls. (1-13) 
Similarly for z, we get 

_ _jcx V 2 z r^ = o i n the bulk 
4 

{d n + jd t )z(r) = at the reflecting walls. (1-14) 

At the absorbing boundaries, we simply impose z = z = 0. These equations are equivalent 
to the zero frequency limit of those in Sec. IB. We observe that the boundary conditions on 
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z and on z are not consistent, which means we cannot find any nonzero configurations z(r) 
that satisfy both these conditions simultaneously. This is due to the fact that the differential 
operator, that appears in ( |1.12| ) between z and z, is not self-adjoint. As noted in KY and 
ML, equations similar to eqs. (|1.13|) and (|1.14p define the right and left eigenfunctions of this 
operator; these eigenfunctions are not complex conjugates of each other. The eigenfunctions 
are not simple to obtain in our geometry, because the x and y dependence does not separate, 
when 7/fl. The analysis of the eigenfunctions in KY and ML ignores the boundary 
conditions at x — 0, L, and is appropriate only for L/W — > 00. 

To obtain a (zero-frequency) diffusion propagator, we use So as the action in a functional 
integral, and define d(r, r') by 



d(r,r') = ^f((z(r)z(v'))) = ^f J D[z,z]z(r)z(r')e s ^/Z . (1.15) 

(Here Zq is the same functional integral without the z(r), z(r) inserted.) Then d(r, r') 
satisfies 

- V 2 d{r, r') = -V 2 d{r, r') = S(r - r') in the bulk; 

(9 n — ■yd t )d(r, r') =0, r at the reflecting walls; 

(& n + r )& t )d(r, r') = 0, r' at the reflecting walls. (1-16) 

The propagator exists and can be shown to satisfy the stated conditions. One notices that 
the propagator is not symmetric with respect to r and r' since the boundary conditions for 
r and r' at the edges differ by a sign. In principle, the propagator can be evaluated by 
expanding in the right and left eigenfunctions, as in KY and ML, however as these are not 
readily available for our geometry, we will just define it by eqns. (|1.16| ) (see also Sec. HID 
below) . 

The action ( |1 . 1 1|) has also appeared in the literature in connection with an open string 



with opposite electric charges attached to the ends, in a uniform magnetic field ||48|| . The 
boundary conditions have also appeared there, along with explicit results for the diffusion 
propagator d in some geometries simpler than ours. We will see later (in Sec. IE and in Sec. 



Ill) that ( |1 . 1 1 ) also arises as the lowest-order part of the nonlinear sigma model action. In 



fact, the full Chalker-Coddington model with phase coherence [43|] is related to the nonlinear 
sigma model [^9| in a manner closely analogous to the relation between the models discussed 
in Sec. IB and here, which are just the linearized versions. In Sec. Ill, we will also discuss 
the expressions for the currents, like those in Sec. IB, from the point of view of the nonlinear 
sigma model action. 



D. Bilocal conductivity tensor and conductance 

In Sec. IA above we have used a classical formulation of the conductivity and worked 
out some of the consequences for the two-probe conductance in high magnetic field. In this 
subsection we introduce a full quantum formulation for the bilocal conductivity tensor and 
the conductance in a disordered phase- coherent system in order to treat both the average 
quantum conductance and its variance. 
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For a quantum conductor with phase coherence, that is, when Li n is larger than the 
sample size, the electron wave function is sensitive to the external field in the entire space. 
Equation ( |1 . 1| ) cannot be used as there is no definition of the chemical potential within the 
sample. We apply standard linear response theory to a finite disordered region (denoted by 
A) with Fermi energy Ep, connected to perfect leads held at fixed voltages, which induce 
a local electric field in the disordered region the detailed form of which is not relevant. 
(For the two-probe case, the sample occupies the region < x < L, < y < W; see 
Figure [I].) Following the treatment of Baranger and Stone |^2| one finds that there is a 
non-local relation between the current response and the applied electric field: 



3i*( 



h J a 



d 2 r' a,Jr,r')Ejr') 



;i.i7) 



where the bilocal conductivity tensor a(r, r') (which has dimensions of inverse length 



squared) at T = can be expressed in terms of a pair of Green's functions [42 



<V(r,r 



where 



and 



Ami 



G + (v,v';E F ) DlD u G-(v',r;E F ) 



Ami 



dE' 



d 



dE' 



-G + (v,r',E') D^D V G + (v',r,E') 



+ G-(v,v',E') dId u -^jG-(v',v,E') 



G±(E) 



1 



E-H±iri' 



X18) 



G{r', r) D G(r, r') = G{r', r)[V - i{e/h)A {r)]G{r, r') 

- G(v,r')[V + i(e/h)A (r)]G(r',r). 

Here H is the Hamiltonian, discussed further in Sec. II, A is the vector potential repre- 
senting the background magnetic field, and we will also use —i%T) = —ihV — eA . In the 
presence of the magnetic field the bilocal conductivity tensor is not entirely a Fermi-energy 
quantity. Even at T = 0, the complete current response function, in the presence of magnetic 
field, contains not only terms involving G + G~ at the Fermi energy, but also terms involving 
G + G + and G~G~ integrated over all energies E up to the Fermi energy. We denote the 
G + G~, G + G + and G~G~ terms as cr+~, cr++ and <J~~ . In disordered systems, products of 
Green's functions that are both retarded or both advanced are generally short-ranged be- 
cause of the amplitude cancellations among different wave fronts (they typically only extend 
over the range of the mean free path), so we can treat cr ++ and a as contact terms: 



a aa (r,r' 



a aa S(r-r'), 



1.19) 



where a = +, — , a aa = J A d 2 r a aa (r,r'). 

In the presence of a magnetic field B, the current given by ( |1.17p does not necessarily 
satisfy V ■ j = even when E is time-independent, unless we also require V x E = 
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(i.e. that the component B of the magnetic field perpendicular to the 2D layer is also 
time-independent). As shown by Baranger and Stone [fEql , under this condition "current 
conservation", V • j = 0, is satisfied, and from it we can derive conditions on the bilocal 
conductivity. Writing E(r') = — V'<f)(r') (where <fi is the electric potential), we have 



j(r) = jj d 2 r> a(r, r')- V' 0(r') - j £ 0, jf _ a(v, r') • dS', 



;i.20) 



where 0j is the (constant) potential in the ith lead, Cj is a surface across the zth lead, and 
the boundary terms at the reflecting walls vanish because the normal component er n „(r, r') 
vanishes for r at the wall (and similarly for v = n, r' at the wall). Thus V • j = implies 
that the following conditions are satisfied: 



V -o-(r,r')- V = 

V- / <t(t, r') ■ dS' = foralH. 

Ja 



(1.21) 



The above identities have been verified in Ref. j|2] for the exact bilocal conductivity tensor. 

Using the second identity in ( |1.21|) , one can transform the linear response equation (|1.17|) 
into a different form. Assuming, without loss of generality, that \x = Ep in all the leads (we 
always view Ep as a position-independent constant), the electric potentials there differ from 
the voltages (electrochemical potentials) only by a constant, and the total current in the zth 
lead can be written as a function of only the voltages in the leads: 



h 



3 1 



'1.221 



where Qij S 3X6 conductance coefficients. The related to the bilocal conductivity 

tensor by 



9ij 



dSi 



a r, r 



(1.23) 



where Si and Sj are cross-sections in the ith and jth lead, and dSi, dSj are differentials 
of outward-pointing normals. (In eq. ( |1.22| ), we used the relation J2j 9ij = 0, which follows 



and implies that a constant V produces no 



from eq. (|1.23|) and the second of eqs. ( |1.21[ 
current in any lead.) For cross-sections Si not intersecting Sj, for i ^ j, the off- Fermi-energy 
terms a ++ (r, r') and a (r, r') can be shown to give zero contribution [42|, and therefore 
can be expressed as a Fermi-energy quantity. It has been shown that is proportional to 
the total transmission coefficient of the scattering states at the Fermi energy from the ith 
lead to the jth lead |4l|j50|j42]1 . In this article we will consider only the simple case of the 



two-probe conductance, in which certain further simplifications are possible. 

Since the total currents at all cross-sections are the same for a two-probe set-up, we can 
average over all cross-sections to get a volume-integral form for the two-probe conductance: 
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d r (T-Jr, r 



1.24) 
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Here we must be careful to include the off-Fermi-energy terms since they are needed to 
preserve the uniformity of the current across each section (from the technical point of view, 
since x and x' can now coincide, we are no longer justified in dropping the a aa terms). 
Despite this disadvantage the volume-integral form of g is often more convenient to use 
in actual calculation than the surface-integral form since volume-averaging still eliminates 
many diagrams which would be non-zero in the surface-integral approach. Thus equations 
( |1.23| ) and ( |1 . 24| ) will serve as the starting points for the evaluation of quantum 



conductance and conductance fluctuations. 

We may gain further insight into the meaning of the current conservation conditions, and 
the relation to the classical case, by use of the self-consistent Born approximation (SCBA) 
results for the disorder average of the bilocal conductivity tensor. As we will show in Sec. II, 
in this approximation (cr^(r, r')) (the single angle brackets will always denote the average 
over the disorder) is of the following form when r and r' are more than a mean-free path 
from the edges: 

<<V(r, r'))scBA = [a°„ 8, v + (a^° + <- ) e ^]5(r - r') 

" 4-[4A + <4VfyH<W - ^V^,]d(r,r'), (1.25) 

XX 

where a xx , cr^°, cr^'°, a xy ° + a xy '° = °% are the SCBA conductivity parameters a**' 
comes from the a ++ and a part of equation ( |1.18| ), and d is the diffusion propagator 
discussed in Sec. IC (with a xy appearing in the boundary conditions). In the zero field limit, 
oil? = al 1 ,; = 0, the above reduces to 

y j, y i 

Mr, r'))scBA = ct° xx [5,J(v - r') - d.&J^r, v>) 

(in which d° is the diffusion propagator for a xy = 0). To our knowledge this basic result first 
appeared in the mesoscopic physics literature in ref. [ 58|| , although it may well have been 



known earlier. It or (|1.25 ) show that the current response to an electric field has a nonlocal 



part, the term containing d° or d in the formulas, due to diffusion. For non-zero magnetic 
field, one finds from eq. Ql.25 ) that 



<9 M (<V(r, r'))scBA = e^Jfd^v - r'). (1.26) 



The divergence of the response current is therefore 



V ■ (j(r)> = jj d\' V • (a(r, r')) SC BA ■ E(r') = £<#°V x E(r). (1.27) 

In the presence of a magnetic field, the current is divergenceless only when V x E = — = 

0; otherwise there is a time dependence in the density, dp/dt = (e 2 / h)a X y ,0 dB jdt. To 
obtain a truly static response, we would have to impose V x E = 0. (There is also an 
edge-current contribution involving a xy , which will be described in the SCBA case in Sec. 
II.) This behavior is typical of the quantum Hall effect, in which a^e^S^r — r') is the only 
part of (cr(r, r')) that is nonzero in the interior of the system on scales larger than the 
localization length £, and a xy is quantized to integer values. It is the local expression of the 
gauge-invariance argument pi 
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The measured experimental quantity is the two-probe conductance. It is straightforward 
to show (see Sec. II), from an equation similar to (|1.25|) , that the two-probe conductance 
within the SCBA, g°, can be written in terms of the diffusion propagator as: 



9° = -^ x j o dy' j Q dy + 7 ^)(^-7at)d(r,r0. (1-28) 

where x ^ x' are arbitrary. We note that although (cr^r, r')) in the bulk depends on cr^y, 
the conductance depends only on the full cr® y , due to additional edge current contributions 
to a which we omitted in ( |1.25| ). These contributions are similar to those discussed in Sec. 
IB. 

It is in fact possible to show that the mean bilocal conductivity and conductance obtained 
in SCBA are identical to those obtained from the "classical" formulas of Sec. IA, if all' is 

I xy 

zero. To obtain the response to an arbitrary electric field E, we write the conditions of Sec. 
IA on the current density, using £ = E + V/i/e, as 

V 2 /i/e = -d^Ep + ie^E v ) (1.29) 

and the boundary conditions 

(d n - 7 5tWe = ~{E n - 1 E t ) (1.30) 

at the reflecting walls, fi/e = Ep/e in the leads. These inhomogeneous equations for /i/e 
are equivalent for 7 = to a 2D electrostatics problem with a mixed Dirichlet- Neumann 
boundary condition |52], and can be solved using a Green's function technique. A slight 
generalization of the same technique works for 7 7^ 0. The required Green's function is 
precisely d(r, r') as defined in Sec. IC, and one finds 

Hir)/e = E F /e - J d 2 r' [(&„ - 7 e^^,)^(r, v')]E v {v'). (1.31) 



Using eqs. fll.lj) , the bilocal conductivity tensor that results is exactly of the form ( |1.25 ), with 
a xy ° = 0' ari( ^ 110 additional edge contributions. Consequently, the two-probe conductances 



g° given by Eqs. (|1.2j ) and (|1.28| ) are the same, for the same values of cr xx and cr° , and 
since this involves only the total cr° , it remains true even if 7^ is included as in 
the SCBA. This implies that the bilocal conductivity in SCBA has just the form which 
follows from a local relation between current and electromotive field, even though there is 
no sensible definition of a local chemical potential in the phase-coherent limit. Hence all the 
conclusions drawn in Sec. IA can also be applied within the SCBA. 



In general, there is also an edge-current contribution, with coefficient cr X y'°, to eq. ( |1.25| ) 



which will be described in Sec. II. There is no reason why both the bulk and edge cr X y 
contributions should not also appear in the so-called classical formulation of Sec. IA, even 
though they were not included in Ref. [j39" |, since terms of this form would presumably still 
be present even if there were inelastic scattering. 

We may also connect the results of Sees. IA and ID with the field theory in Sec. IC (again 
for a™' = 0). In Sees. IB and IC, the external electric field was zero. If we replace d^z, 
d^z by d^z — 2iA^, d^z + 2iA^ in the action S , eq. (this A should not be confused 
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with A or any other "physical" vector potential) and then define j M (r) = 6So[A]/6A li (r) 
to be the current, we find = ia^ u (d u z + 2iA u )/2, and the conditions for an extremum of 
the action (equations of motion) are d^j^r) = in the interior and j n = at the reflecting 
walls. Thus, these equations have the same form as the local classical conductivity equations, 
with iz/2, —A^ in place of (/z — Ep)/e, E^, and since, for a quadratic action like So, the 
linear response 5((j fJ ,))o/5A u to the perturbation A can be equivalently obtained either from 
calculation of the correlation function (similar to that in the definition of d, eq. ( |1.15| )), or 
using the equations of motion as above, we get the same bilocal conductivity in this approach 
also. This calculation is done in detail in Sees. Ill and IV, so we refrain from giving further 
details here. Sec. IIIC also includes the corresponding "classical" equations for the case 
when erg' ^ 0. 

The SCBA for the average conductance is the leading approximation in an expansion 
in powers of l/a® x , and it is not really surprising that this leading approximation behaves 
identically to the classical case. When evaluating conductance fluctuations, or weak local- 
ization corrections, one must consider higher orders in l/cx^,. In such calculations the tilted 
boundary condition is modified further. In the framework of diagrammatic perturbation 
theory, this can be alternatively viewed as the appearance of additional boundary vertices 
describing interference effects. These vertices are more easily obtained and evaluated in the 
nonlinear cx-model approach to which we now turn. 



E. Nonlinear cr-model approach 

The approach of the previous subsection, in which the self-consistent Born approxima- 
tion is the leading contribution to conductivity and conductance, can be developed as a 
diagrammatic expansion (see Sec. IIA). However this approach becomes cumbersome in 
higher orders because all diagrams contain vertices which need to be evaluated in terms 
of the average single particle Green's functions, and are dressed with disorder lines in all 
possible ways. However, when these vertices, which describe interference or "interactions" 
of diffusing modes, are calculated at small wavevectors, they are all found to be related to 
the same quantities cr° x and a xy . These complex and often redundant calculations can be 
avoided by using the nonlinear cr-model (NLcrM) representation of the problem. 

The NLcrM approach starts by considering only Green's functions at the Fermi energy, 
which means that the non-Fermi-energy parts of the conductivity described in the previous 
subsection cannot necessarily be obtained, though the conductance and its fluctuations 
can. After introducing replicas, the disorder is integrated out, followed by the variables 
representing the electrons propagating at the Fermi energy. After a Hubbard- Stratonovich 
decoupling, and neglecting modes that have no long-range effects, one is led to the action 
^8| , |30[| (more details will be given in Sec. Ill) 



S = f dVtr^Q^Q] - ^ / dVtr^Q^Q&Q], (1.32) 

o J A oJA 

where Q is a 2n x 2n Hermitian matrix obeying Q 2 = I 2n , (l2n is the identity matrix); Q has 
n eigenvalues equal to +1, n equal to —1. Using a parametrization given explicitly in Sec. 
Ill, it can be shown that for small fluctuations about the maximum action configuration 
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where Q is diagonal, S reduces at quadratic order to n 2 copies of the earlier action So in 



Before discussing perturbation theory for this action, we wish to mention some general 
issues. The second term in S is the so-called topological term. It is the only possible term 
that can be added to the first term in two dimensions that is consistent with the symmetries 
of the problem and contains only two gradients (higher derivatives would be irrelevant at 
long length scales). On a compact, oriented manifold without boundary, such as a sphere or 
a torus (i.e. periodic boundary conditions), this term (with a factor of 2ma^. y removed) is a 
topological invariant, which takes integral values. This follows from the fact that the term is 
a total derivative, and the absence of boundaries. Consequently, only the value of a® y modulo 



1 is important. Moreover, because the term is a topological invariant, a does not appear 



xy 

in perturbation theory at all, but only in non-perturbative effects involving configurations 
( "instantons" ) for which the topological term is nonzero. However, in the integer quantum 
Hall effect, we expect to obtain plateaus at integral a xy and transitions between them at a xy 
half-odd- integral (a xx , a xy denote the renormalized, large-scale, parameters, as opposed to 
the bare values a xx and a xy at the cutoff scale I). Because of the periodicity in a xy , the NLcrM 
predicts that all these plateaus and transitions will have identical universal properties. But 
by the same token, it is also unable to predict the integral part of a xy that would be observed 
in a measurement; this information appears to be lost in going to the NLcrM. 

The apparent paradox is resolved on examining the action S for a system with reflecting 
boundaries. The "topological" term is a total derivative that can be rewritten as a boundary 
term, just as for the action So above. The boundary term is not a topological invariant, so 
it can affect perturbation theory, and since it takes arbitrary real values, the magnitude of 
a xy is important, and not just its value modulo integers. Thus when the boundaries of the 
system are correctly taken into account, the value of (J X y CcLll be obtained within the NLcrM 
formulation. Some additional remarks about this point are made in the Conclusion. 

Since the leading order part of the action is the same as So, the propagator for small 
fluctuations in Q is the same as the propagator d discussed earlier, and depends on o\ y 
through the boundary condition. In the work of KY and ML, this modified propagator was 
the only effect included, and was just inserted into the Xiong-Stone results for conductance 
fluctuations. However, the NLcrM shows that the boundary term also contributes at higher 
order, producing new vertices for "interactions" between diffusons, which are boundary 
interactions with coefficient a xy , and which contribute to the fluctuations at leading order. 
These terms must be present in order to maintain the full U(2n) symmetry of the NLcrM, 
which essentially corresponds to preserving the continuity equation for the current. We 
have also obtained them in the diagrammatic approach, but only with much additional 
effort. (Pruisken [30[] also discussed the relation of the topological term to edge states, but 
appears to infer an incorrect boundary condition. His boundary condition is very useful in 
instanton calculations p9|-pT|1 but does not correctly represent the edge effects, unlike the 
conditions to be discussed in this paper.) In this paper, we evaluate the effects of these terms 
to leading order in \ja\ x , and, to simplify the calculations, also to leading nontrivial order 
in 7 = cr xy /a xx . As well as calculating the mean and variance of the two-probe conductance, 
we show how the expression for the mean bilocal conductivity tensor in the SCBA can 
be recovered within the NLcrM, including the non-Fermi-energy parts, by modifying the 
coupling of the NLcrM to the external field, and we discuss the resulting form of current 
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conservation conditions to all orders in perturbation theory. 

One may wonder if the boundary conditions for a system with boundary invalidate the 
conclusions of the analysis of Pruisken and co-workers |29H3Tf, who studied effects of instan- 
tons in a system without a boundary. Strictly speaking, in a finite system with boundaries, 
there are no well-defined topologically distinct sectors. However, small instantons which 
are well localized inside the system boundary, so that Q approaches a constant outside the 
instanton core and satisfies the boundary conditions at the edges, while probably not exactly 
local minima of the action, are still nearly so when their size goes to zero compared with 
the system dimensions, and in this limit their topological charge will still be an integer, and 
they will make the same contribution to the action as they did for the other b.c.'s. Thus, 
the effects on renormalization group flows for a® x and a® y , obtained in the interior of the 
system, should be unchanged. These effects will not be considered further in this paper, 
which emphasizes perturbation theory. 

In the one-dimensional case of the NLcrM, no other term can be added to the basic 
gradient-squared term, and there are of course no sides on which the boundary condition 
could be modified. Therefore, the conductance fluctuations in the unitary ensemble should 
be universal, and must be recovered in the quasi- ID limit of the 2D system in a magnetic 
field which we are considering; this constitutes a strong check on the 2D calculations. We 
will show that in this limit the only effect of a xy is to modify the ID conductivity, which 
is known to cancel from the conductance fluctuations to leave a universal number. KY 
and ML claimed that a xy does affect the quasi- ID limit; the present argument shows that 
their results must be incorrect. In the Conclusion, we will briefly mention the situation for 
dimensions higher than two. 



II. DIAGRAMMATIC EXPANSION FOR (cr^(r,r )) 
In this section, we evaluate the mean bilocal conductivity tensor for the short-ranged 



potential model using the diagrammatic impurity- averaging technique [53]. We will first 
review the self-consistent Born approximation (SCBA), which is the leading order of the 
perturbation expansion, and establish basic parameters such as the mean free path I, and 
the bare conductivities a® x and a° y . A gradient expansion is used to treat the current vertex. 
Within the SCBA, (cr^r, r')) has a contact term as well as a long-ranged term which can be 
expressed in terms of the diffusion propagator. In the bulk, a!£ appears in the long-ranged 
term. Along the reflecting boundary, the edge currents give rise to 5-function contributions 
proportional to cr^'°. As a result, it is a° y = a^' y + a xy '° that appears in the boundary 
condition and the two-probe conductance. We will also check that current conservation is 
respected within the SCBA. 



A. The model, edge current, and ideal leads 



An electron in a system with edges, in a random potential and subject to a perpendicular 
magnetic field is described by the following Hamiltonian (we neglect spin throughout this 
paper): 
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H = H + V(r), H = ^-(-thD) 2 + U(y), (2.1) 

where V^(r) is the random potential and it is confined in the region of < x < L, and U(y) 
is the confinement potential (see Figure [I]). The uniform magnetic field is in the z-direction, 
and we choose the gauge A = —Byx. For simplicity, let us assume the confinement U(y) 
to be the hard wall potential with U(y) = for < y < W, and U(y) — oo for y < 
and y > W. The infinite potential barrier requires the wave function to vanish at y = 0, 
W. The system is infinite in the ^-direction, but the disorder is present only in the region 
< x < L. For the random potential, we will take the simplest among all short-ranged 
models, which has the statistics of the Gaussian white noise with zero mean (angle brackets 
denote the disorder average): 

(V(r)> = 
(V(r)V(r')) = uS(r-r'), 
(V(r l )V(r 2 )---V{r n )) connccted = 0, n > 2, (2.2) 

where u describes the degree of disorder. 

In the absence of the random potential, the unperturbed wave functions can be found 
by separating the variables [10]. The wave functions ip{r) are labeled by the wavevector k 



in the x-direction and by N in the transverse direction; N turns out to be the Landau level 
index. We have ^ 

^TVfc(r) = -j=e ikx (j) Ntk {y), 
V L 

and (f>N,k{y) satisfies: 

{% 2 [{-id y f + (y- l 2 B k)H B A ]/2m e + U{y)}<t> N , k {y) = E Nik (f> Nik (y), (2.3) 

where l 2 B = H/eB. Without the confinement potential, the Hamiltonian is simply that of 
a harmonic oscillator, with the harmonic potential centering at y k = kl 2 B . We have [50 



4>N,k{y) = XN(y - l 2 B k), E Ntk = E N = (N + l/2)huj c for W > y k > 0, where xn is the 
iV-th wavefunction of the harmonic oscillator, and = 0, 1, 2, .... The wave functions 
spread an extent R c = \/2N + 11b around y k . We can see that for wave functions which 
center at a distance more than R c away from the walls, the presence of the walls is hardly 
felt, but for those which reside within a distance R c from the walls, their eigenenergies are 
raised above En, because the wave functions are forced to zero at the boundaries and thus 
made to oscillate more rapidly near the walls. Only the states within R c of the edges have 
non-zero group velocity along the walls, i.e., the expectation value of the velocity operator 
(N, k\v x \N,k) ^ 0. From now on, R c plays the role of the short length scale of the problem 
and we treat the edges as having zero width. We will later show that, in this sense, the 
inhomogeneity at the edges gives rise to 5-function contributions to the current (see also 



Although the above description of H is very convenient for finding an explicit solution 
for the energy eigenf unctions when the system is infinitely extended in the x direction, 
it does not provide a convenient description of ideal leads in the presence of a non-zero 
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magnetic field. For our purposes ideal leads should be perfect absorbers of all incident 
current, i.e. they should behave as if they have essentially infinite conductivity compared to 
the sample. The problem is that in the leads, where x < or x > L, the states at the Fermi 
energy generically consist of a certain number of edge channels moving in each direction, 
and this number is equal to the number of Landau levels below the Fermi energy in the 
bulk, which is N + 1 when the Fermi energy lies above the iVth Landau level. Thus in the 
relevant sense, the leads have zero bulk conductivity and only absorb and inject current at 
the corners. There are several ways in which we could modify our model to remedy this 
problem. One way would be to let the magnetic field drop to zero at the ends of the sample, 
so the leads are like the usual 2D metallic leads in zero magnetic field. Another way, which 
corresponds roughly to the ohmic metallic contacts used in real experimental systems, would 
be to "thicken" the system outside the sample so that it becomes three dimensional, thus 
increasing its conductance. There is however a third way, which most closely conforms to 
the perfect leads used in the network model (see Fig. ^) and has some convenient properties. 
The links of the network can be viewed as edge channels, and outside the sample there are 
many of them, running parallel, alternately right- and left-moving, without backscattering. 
A similar setup can be produced in a 2D Hamiltonian model with a uniform magnetic field, 
by replacing U(y), only in the leads, by a potential Ui(y), which = +00 for y > W or y < 0, 
and has a sinusoidal form in < y < W. If the Fermi energy lies between the maxima and 
minima of Ui, there will be many "internal" edge channels at the Fermi energy, consisting 
alternately of N + 1 channels moving in one direction and N + 1 moving in the other. For an 
infinitely-long, translationally-invariant system there will be no backscattering among these 
modes. In effect, we have many narrow leads in parallel, all connected to a single reservoir at 
—00 and to another at +00. Then the number of right-moving channels can be proportional 
to the width of the system, or arbitrarily large, and the current can be injected uniformly 
across the end of the sample. From now on, it is this model that we will implicitly use. 



B. The self-consistent Born approximation (SCBA) 

In order to calculate (a M ^(r, r')), we first need to evaluate the disorder-averaged single- 
particle Green's function and two-particle Green's function, (G) and (GG), which can be 
expanded in power series in l/{kpl). The SCBA takes into account all the non-crossing 
diagrams (Figure and it has been shown to be the leading contribution in l/(kpl) f27fl . 

Within the SCBA, the single-particle Green's function (see Figure 0(a)) satisfies 

[E - ifo(r) - S(r)]G(r, r'; E) = 8{v - r'), (2.4) 

where the self-energy in turn depends on G (we use the over-bar to denote SCBA Green's 
function): 

£(r,£) = uG(r,r; E). 

The above equation can be solved analytically in the low field limit (uj c t <C 1) and the high 
field limit (uj c t ^> 1) |24]]. In the intermediate field range, it can be solved numerically. 
For our purpose, we do not need the explicit solutions. We use the SCBA Green's function 
to define the effective scattering rate 1/r and the effective mean free path / at the Fermi 
energy: 
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1 _ £ + - E- / 1/ro, ^ c r «l , , 

r ^ 1 2rsin^ F )/^, ^ c r >l ' 1 - 0) 



where T = ^Ju/(2ttI%) ~ %Juj c /tq is the width of broadened Landau levels, 6(E) 
< 6(E) < 7T and 



cos 1 



E-E N 
2T 



I 2 = u f d 2 r x (x - x 1 fG + (r,r 1 ,E F )G-(r 1 ,v;E F ) = $ << j . ( 2 .6) 
j y k c} uj c t ^> i 

Here Iq is the mean free path in zero magnetic field, 1$ = hkpTo/m e . 

Within the SCBA, the two-particle Green's function S + ~(r,r') = uG + (r, r')G~(r', r) 
amounts to adding up all the ladder diagrams (figure 0(b)). The sum can be written in the 
form of an integral equation: 

S+-(r,r', £,£') = 

uG + (r, r', E)G~(r', r, E') + J d 2 r x uG + (r, r 1; E)G~(v x , r, E')S + ~(r 1 , r', E, E'). (2.7) 

We make use of the fact that G (r, r')G (r',r) is short-ranged and expand S' + ~(r 1 ,r / ) in 
the vicinity of r t = r. We get, for E, E' close to E F , 

J d 2 r 1 uG + (r,r 1 )G-(r 1 ,r)S+-(r 1 ,r') = 

[C (r; E, Ef) + C a (r; E, E')-V+ ~C 2 (r; E, E')V 2 + ■ ■ -]S + -(r, r') (2.8) 



where 



C (r; E, E')=u[ d 2 n G + (r, n, E)G (r x , r; E') ~ 1 + i(E - 
Ci(r;S,^) =« | rf 2 n (r-rOG+^n;^)^^,^ £") ~ 0, 
C 2 (r; E')=u J d 2 n (x - £i) 2 G + (r, n, £)G~(r x , r ! E ') = ? ■ 
It follows that S + ~(r, r', E — E') satisfies the diffusion equation: 

[-VtV 2 - i(E - E')r/h}S + -(r, r; E - E') ~ J(r - r) (2.9) 

where D is the diffusion constant: T>(E) = l 2 /(2r). One can see that at E = E' = E F , 
which is all that will be required in this paper, is proportional to the dimensionless diffusion 
propagator d(r, r') we defined in equation ( |1.16| ): 

^ + -(r,r';0) = d(r,r'). 

We postpone derivation of the boundary conditions on d until after we have discussed the 
conductivity tensor. The ladder sum for S ++ = uG + G + and S = uG~G~ can also be 
carried out in similar fashion. It is easy to see that S ++ (r,r') and S (r, r') are generally 
short-ranged. 
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Now we are ready to treat the mean bilocal conductivity tensor. Within the SCBA, 
(er Mi/ (r, r')) has two contributions: the simple bubble diagram and the sum of the ladder 
series (Figure 0(c)): 

(<V(r, r'))scBA = <v(r, r') bubb i e + cv(r, r') 

ladder • (2-10) 

The bubble diagram has the range of the mean free path I and we treat it as a ^-function: 

a^(r, r') bubblc = orJ„5(r - r'). 
cr°„ are the SCBA conductivity parameters, which are essentially constant inside the sample: 

a °uu = f dVcr^(r, r ) bubWe = 7777 f d 2 r [ d 
J A LW J A J A 

It follows from the definitions (c.f. ( |1.18| )) WW>^ that 



2 r'cr^(r, r') bubb i c . 



(2-11) 



where AG = G 



-Tr 

LW 

LW 

h 2 

2LW 
G~ 



v x G^(E F )v x G (E F ) 



d 



dE^Tr 

Tr {v x AG(E F )v x AG(E F )} 



Tr 



v T G u. 



d_ 

dE 



G~ 



(2.12) 



-ihD/m e , and Tr denotes the trace of the matrix product, in 



which G (r, r') are viewed as r, r' matrix elements. For a 



xyi 



JO 
a xy 
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1,0 

xy 



1,0 1 11,0 
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h 2 



LW 



:Tr 



v x G + (E F )v y G (E F ) 



a 



11,0 

xy 



LW 



E F 



dE' I Tr 



v x 



_d_ 

~dE' 



G+ \v„G 



+ Tr 



v x G v % 



_d_ 

dE 



■G- 



(2-13) 



The above expression for a xy can be put in the same Fermi-energy form as in Ref. |55 



and a® have the following limiting behavior |]24]p8| , |30|| : 



a° xx = hV(E F )p(E F ) 



hn e T /m e , u c t < 1 

(2N + 1)-k~ 1 sin 2 6(E F ), co c t > 1 



a 



xy 



(T° X LU C T , UJ C T < 1 
N + V, UJ C T > 1, 



(2-14) 



where p(E) is the SCBA local density of states, n e is the electron density, N is the highest 
Landau level index, N = 0, 1, . . . , and v = 1 — 6/tc is the filling fraction of the highest 
Landau level, which is well-defined since there is vanishing local density of states in the bulk 
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in between the Landau levels in the regime oj c tq ^> 1, within SCBA. Often, only the peak 
value (2iV + 1)/tt of <r® x is quoted in the literature; we emphasize that a xx has oscillations 
and goes to zero when the Fermi energy lies in one of the gaps in the bulk density of states. 

Now we treat the ladder diagrams. Since S ++ , S are short-ranged, a ++ and a are 
also short-ranged. One can further show that the ladder series for and a~~ do not give 
additional contributions to the 5-function term. The long-ranged term of <7^(r, r') comes 
from 

<7j7(r,r')iaddcr = h 2 u J d\ x J d 2 r 2 J+"(r, r 1 )S' +_ (r 1 , r 2 ) J~ + (r', r 2 ). 
The current vertex J(r;ri) is short-ranged: 



r b (r, ri ) 



<n 5 6 (r X) r)D5 B (r ) r 1 ), 



2m 



(2.15) 



where a, b = +, — , and D was defined in Sec. IC. We can carry out the integral by expanding 
«S ,+ ~(ri, r') in the neighborhood of r: 



J d 2 n J(r, ri )5 + -( ri , r 2 ) = [J (r) + J x (r) • V + - • •] S + -(r, r 2 ). 



The first term is 



-ih 



2m P 



1 



d 2 r x G (r 1 ,r)DG + (r,r 1 ) 



-(r|v(G + -G-)|r). 



(r|v(G + -G )|r) as 2vrz^ 5(E F - E a )^ a (r)v^*(r) 



(2.16) 



In the last step, we have used the identity G + G = (G + — G )/ (S + — £~). Rewriting 



(2.17) 



it is easy to recognize that Jg~(r) is proportional to the impurity-averaged equilibrium 
current density at position r and at the Fermi energy Ep- In the bulk, this should be zero 
by isotropy, however, at the edges isotropy is broken, Jq ~ parallel to the boundary is not 
zero. Hence we can write 



J 0,: 



2ni 



e(£- 



dI e (E) 



dE 



S(y -W)- 



dI e {E) 



dE 



Kv) 



where I e (Ep) is the total edge current. It has been shown in Ref. [30] that 



Oh 

dE 



dM 



E F 



dE 



11,0 

xy 



Ep 



h 



where M(E) is the total magnetization. 

The coefficients for the first-derivative terms are: 



[2.18) 



[2.19) 
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n 



- / " / LW 
1 



— f Tr[vpG + r v G } - Tr[r u v fl G + G } + ^-Ti[G + G ]sA, 



LW 



Tr[v„G r v G ] - Tr[rVG G ] + ^Tr[G G ]5 



Expressing the coefficients in terms of the SCBA conductivities and the mean free path, we 
have: 



u fo J 1, xx J t, xx — a xxl /2j 

M ^ J l, xx J l, xy = a xyl 1 2- 

Putting all the pieces together, we obtain for the full bilocal large-scale conductivity 
tensor within the SCBA, including the edge effects: 



(o>(r,r'))scBA 



11,0 

xy I^V 



5{r 



a'' 



+ + - W) - %)) 



x 



°IA - - -W)- %'))! d(r, r'). (2.20) 



Within the SCBA, the conductivity tensor obeys (<7 yi/ (r, r')) = 0, for r at the reflecting edges 
and r sufficiently far from r', as an exact relation before the long- wavelength approximation 
is made. Together with ( |2.20 ), this implies the large-scale boundary conditions ( |1.16 ) on 
the diffusion propagator d; see Appendix A. Because of the edge current, this is very similar 
to the discussion of the conservation of the j' current at the edge in Sec. IA (where oj£'° was 
-1). 



C. The SCBA for the two-probe conductance 



We are finally ready to express the two-probe conductance in terms of the surface inte- 
gral of the diffusion propagator. For the integrated currents at any cross-sections, the two 
opposite edge currents at the boundaries can be transformed to an additional bulk derivative 
in the transverse direction: 



ir 



dy[5(y-W)-5(y)]d(r,r') 



dy d y d(r, r) 



therefore, the second and third terms in the square brackets of ( p. 20 ) combine in this case 
to give o~® y . We get the form for the SCBA conductance g° (|1.28| ) expressed in terms of an 
integral over two cross-sections. This can be further transformed into the following volume- 
integral form: 



9° = ^ 2 J A d 2 r l/ 2 r'al[5(r-r')-(d x + 7 d y )(d f x -jd' y )d(r,r')). 



(2.21) 



In the next section, we will see that the above expression for conductance can also be 
obtained using the NLaM formalism which we develop below. 
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III. FIELD-THEORETICAL APPROACH 



In this section, we first set up (in Sec. IIIA) a generating function from which we can 
obtain any expressions involving the Green's functions at the Fermi energy. As we have seen 



in the introduction |[42|| , this allows calculation of the mean and the variance of the two- 
probe conductance. In Sec. IIIB, we discuss the coupling to a U(2n) vector potential, which 
can be used to generate the NLcrM conductivity. For a certain form of coupling, we can in 
fact recover, at lowest order in l/u xx , the SCBA form of the bilocal conductivity tensor. We 
discuss the physical meaning of the various terms in relation to current conservation. The 
results show that the variance of the two-probe conductance can be calculated within the 
NLaM. We then turn to the perturbation expansion itself. 



A. The nonlinear a model 



In setting up the partition function (or generating function for average Green's functions), 
we will use the replica method to perform the average over disorder. We define 

Z = J D[V]P[V} J D[ip,ip]exp J d 2 rip(r)[E - H + ir]A]ip(r), (3.1) 

where, as in Sec. II, H = H + V(r), = ...<£>"... {i = 1, 2, . . . , n, a = +, — ) is a 
2n-component vector of complex Grassmann numbers, rj = + , and 




(J n is the n x n identity matrix). Here we remark that the choice of Grassmann number 
(anticommuting) fields, which leads to the symmetry group being U(2n), as will be discussed 
below, is not essential. If one uses bosonic (commuting) fields, the symmetry is the noncom- 
pact group U(n,n). This is usually not used in the Hall effect because there is no topological 
term in this case for n > integer, and so the a xy dependence found from non-perturbative 
instanton effects is not seen JJT]. However, there is nonetheless a boundary term of the 
same structure as in the U(2n) case, and for a system with boundaries this has effects even 
in perturbation theory, and these will be the same in the n — > limit for either choice of 
symmetry. We will continue to work with the choice that leads to the compact symmetry. 
The random potential V"(r) has the Gaussian distribution 

P[V}oce-I d2ry2 ^ 2u , 

as in Sec. II. For rj = 0, the action has global (r-independent) U(2n) symmetry, which acts 
on ip as (p(r) — > Ulp(t), where U is an element of U(2n). For r] > 0, the symmetry is broken 
to U(n)xU(n). In discussing the conducting properties, it is useful to introduce a source 
term that will generate current correlations. This is done by introducing the vector potential 
A(r), where A is a 2n x 2n hermitian-matrix- valued vector field (not to be confused with the 
vector field A associated with the constant magnetic field B = V x A ). It is introduced 
into Z by replacing the covariant derivative D (viewed as multiplied by I 2n ) by D — iA. 
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The generating functional Z[A] then has gauge invariance, since under a local U(2n) gauge 
transformation ip(r) — > Uip{r), — > UA^U^ 1 + Ud^U^ 1 , the action S[A] is invariant, and 
so is the integration measure D[ip,Tp\. Performing the functional integral over ip and Tp, we 
get 

Z[A] = J D[V]P[V]exptrTr\og[E - H -i5A-h 2 (A 2 + i~D ■ A + iA • D)/2m e ], (3.2) 

where Tr denotes the trace over functions in real space as before, and tr is the trace over 
the 2n replicas. 

From the above partition function, we cannot get the exact bilocal conductivity tensor 
cr^r, r') as given in ( [1.18] ), since we cannot generate the below- Fermi-energy contributions. 
However, on length scales greater than the mean free path, the latter simplify and can be re- 
expressed as Fermi-energy terms, as will be shown below, and these can be reproduced from 
our partition function. An exception to this is the cr X y° term in the bulk, which therefore 
describes true non-Fermi-energy physics. We can however get the expression for the two- 



probe conductance, which can be written in terms of Fermi-energy quantities alone 
Let us assume that the source field A is independent of y, then using S/SA(x) to denote a 
functional derivative for a y-independent variation, we can show that 



S 2 Z[A 

— hm 



J D[V]P[V] JdyJ dy' j>TH°g(B-ff-^A) 



(3.3) 

The 5-function term is related to cr+ a + (r, r') and a xx (r, r') (see Sec. ID). We have argued in 
Sec. ID that cr x %(r, r') (a = +, — ) in disordered systems are short-ranged and can be treated 
as contact terms cr™5(r — r'). One can show using the commutation relations v x = i[H,x]/h 
and [v x ,x] = h/im e P2| , |55| , |32| that 



o% = -^Tr [G a {E F )v x G a {E F )v x ] = ^G a (r,v). 



Therefore, the expression in the square bracket in ( |3.3| ) gives an approximate version of 
the unaveraged a xx (r,r'), eq. (|1.18|) , valid on scales greater than the mean free path I. 



The ^-function term drops out for two cross-sections far apart. Taking the limit n — > 0, 
in which case e trTrl °g( E -- ff - ir ? A ) _ » i ; equation ( |3.3| ) gives the disorder-averaged two-probe 
conductance: 

lg(x,x')) = - lim lim - - , _ ^ 5!"^"-+ , a • (3-4) 
n^A^8Al n {x)8A x ^{x') 

Similarly, the second moment of the conductance can be obtained by applying four deriva- 
tives to the partition function: 

(gfx-i, x' 1 )g(x 2 , x' 2 )) = lim lim - ,_ - — , f- ."f [^"! _ , — , r ._, , . « ■ (3.5) 
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These expressions should be independent of xi, . . . , x' 2 . This will be discussed in the next 
subsection. 

An effective NLcrM action can be derived using the Hubbard-Stratonovich transformation 
28| , |30|| ; a mean-field approximation then corresponds to the SCBA. The fluctuations about 



this mean-field theory produce all of the higher-order effects. An effective action for the 
long-range effects can then be derived; we get, retaining only terms with no more than two 
derivatives, and omitting the gauge field A until the next subsection, 

S = J d 2 r {-^tr^Qc^Q] - ^tx^Qd^Q] + r?tr[QA]} , (3.6) 
Z = ( D[Q]e s . (3.7) 



The field Q is a 2n x 2n Hermitian matrix obeying (at each r) Q 2 = J 2n , of which n 
eigenvalues are equal to +1, n to —1. This means that, at each r, Q takes values in the 
coset manifold U(2n)/XJ(n)xXJ(n). The action S has the same symmetry as the original 
action, Q(r) — > UQ^U^ 1 , but Q is invariant under the diagonal U(l) subgroup of U(2n). 
The remaining SU(2n) symmetry is again broken to S(U(n) xU(n)) for non-zero rj. The 
parameters a® x and a® are bare conductivity parameters, like those resulting from the 
SCBA, (but may differ by finite renormalizations, corresponding to short-range effects that 
are not included again by the NLcrM) and describe the response of the system at the scale of 
the short distance cutoff, which is of order the mean free path I. The measure D[Q] = Y\ r dQ 
is the product over points r inside the sample of the unique SU(2n)-invariant measures dQ 
on the space U(2n)/U(n)xU(n) for each point r. At the ends x = 0, L, we insist that Q = A 
to represent the absorbing boundary condition. 



B. Gauge invariance, current conservation, and boundary condition 

In this subsection, we discuss the way in which the gauge potential A enters the NLcrM 
action, and the related questions of current conservation, the equation of motion, and the 
tilted boundary condition. We begin by requiring that the action be gauge invariant. In 
Sec. IIIC, we will modify it to a non-gauge-invariant form to bring the conductivities into 
line with those discussed in the previous sections. 

In view of the gauge invariance of the generating functional Z[A], the action 5 [A], 
including A, should also be invariant (when rj — 0) under the local gauge transformation 
Q -> U{r)Q(r)U- l {r), -> UA^U' 1 - iUd^U' 1 (it is assumed that U respects Q = A 
at the ends; the invariance of the functional integration measure implies that invariance of 
the action ensures invariance of Z[A]). The simplest way to introduce the external source 
field A into the NLcrM is to replace the partial derivative d^Q by the covariant derivative 
Dfj,Q = d^Q + i[An, Q] everywhere in S. This leads to a manifestly gauge-invariant action, 
which is given below as the a® x and a 1 ^ terms in eq. ( |3.8|) . However, this is not the only 
way. The second way is less obvious and will follow a brief digression. 

To obtain the other gauge-invariant coupling to A, let us first point out that the topo- 
logical term, without any A-dependence, is a total derivative, which is a function of the 
values of Q on the boundary and the homotopy class of Q in the interior for given boundary 
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values, but not of the detailed form of Q in the interior. That is, it is possible to change Q 
in the interior, leaving its boundary values fixed, in such a way that the topological term 
changes. Now any change in Q can be viewed as the result of a gauge transformation, and 
a transformation U that leaves the boundary values of Q unchanged must reduce at each 
point on the boundary to some element in a certain S(U(n)xU(n)) subgroup, determined by 
Q at that point. Such a gauge transformation is characterized by an integer- valued winding 
number q, say, which describes the winding of the gauge transformation U at the edge, and 
it changes the topological term by 2nia xy q (for a xy equal to an integer, this has no effect 
on the exponential of the action). In particular, there is a continuously-connected class of 
transformations for which q = 0, so that there are continuously-connected classes of configu- 
rations Q with given boundary values throughout which the topological term takes the same 
value. Further discussion of the topological issues, related to edge states and quantization, 
is contained in Appendix C. 

Now, because the topological term is (apart from the topological effects just discussed) a 
function of only the boundary values of Q, this suggests that we can attempt to compensate 
for a gauge transformation by including a coupling to A on only the reflecting (hard) walls. 
The form of this coupling can be easily obtained, and on including both forms of coupling 
with different coefficients a x ' y , <J xy °, whose sum is a xy = a x ' y + (ji y °, in the respective forms 
of the topological term, we obtain the action which is gauge-invariant when 77 = 0, 

S[A] = J d\ ^- l -alMD,QD,Q) - l -o^tv{^QD,QD u Q) 

+ \°%* J dx [ti(A x (x, W)Q(x, W)) - tr{A x {x, 0)Q(x, 0))] . (3.8) 

The justification for identifying this split of a xy into two pieces cx^ , o"^' , with that found 
in the previous sections, as implied by the choice of notation, will be given below. We note 
that the action is independent of the trace of A. To see this, we may split A into the sum 
of a traceless part and the trace multiplied by J2 n /2n. The latter part is the gauge potential 
corresponding to the diagonal U(l) subgroup generated by lm, and it does not contribute to 
D^Q (because [hnyQ] = 0), or to the cr^'° edge coupling (because tiQ = 0). Hence, S[A] 
is independent of it, but it may be left in for convenience. The action S[A] is easily verified 
to be invariant under any gauge transformation, including the topologically-nontrivial ones 
that leave Q on the boundary unchanged. 

Although we have not given a derivation of our gauge- invariant action S[A], eq. Q3.8| ), 
from the gauge-invariant generating functional (|3.2|) for averages of products of Green's 
functions at the Fermi energy, it is not very difficult to extend the existing derivations (see, 
e.g., |50| ) to include A, and obtain eq. fl3.8|). It is almost self-evident that this will be 
obtained, by comparing (i) the diagrams for the response to A at the Fermi energy that are 
obtained from eq. ( |3.2| ) with (ii) those studied in detail for the mean bilocal conductivity 
tensor within the SCBA in Sec. II, and (iii) those obtained in the perturbation theory for 
the NLuM constructed below. 

We now use the gauge invariance of S in the 77 — ^ limit to derive some current- 
conservation relations for the NLcrM; only infinitesimal, topologically trivial gauge transfor- 
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mations are needed for this. The tilted boundary condition is one consequence of current 
conservation. We will recover the expression for the bilocal conductivity in the SCBA as 
the leading-order term in an expansion in l/a^. x . We will also show that the conductance 



we obtained using equation (|3.4j) is independent of the positions of the cross sections so that 
the volume-integral form for the conductance can be used. 

We begin by considering the equations of motion that follow from the action S[A]. As 
their name implies, these are the equations of motion that are obtained if S[A] is used 
as the action of a classical nonlinear field theory (in one space and one imaginary time 
dimension). The canonical way to obtain these equations is to seek an extremum of S, such 
that 5S/5Q = 0, where the variation, which is the usual one that varies both the x and y 
dependence of Q, respects the restrictions Q 2 = I, Q = Q\ and the boundary condition 
Q = A at the open ends (we note that this is imposed on all configurations in the functional 
integral). The resulting equations also serve as operator relations in the quantum field theory 
that we take to be defined by the functional integral 



Z[A] = J D[Q] e s ^. 



In the functional integral language, the equations of motion become identities among corre- 
lation functions. They are obtained in general by the following argument: Consider a small 
change in Q, Q — > Q' = Q + 5Q, as a change of variable in the functional integral. Since Q 
is integrated over, such a change can have no effect on Z[A]. On the other hand, it changes 
S, and provided the change is such that the Jacobian resulting from the change in measure 
is 1, we obtain the identity 

which is the equation of motion. 

A variation of Q that respects its form can be parametrized as Q' = UQU^ 1 , where 
U = exp iR and R(r) is a 2n x 2n hermitian matrix function of r, and thus is a gauge transfor- 
mation, which leaves the integration measure unchanged; however, we vary Q while leaving 
A fixed. If we view S[A] as a functional of Q as well as of A, thus writing S[A] = S{Q, A} 
(the curly brackets are used to avoid confusion of the two arguments of the functional with a 
commutator), then the equations of motion are equivalent to 5S{Q',A}/5R = 0, evaluated 
at R = 0. This may be re-expressed by making use of the gauge invariance of S{Q, A}. 
Gauge invariance tells us that if Q' = UQU" 1 , A' = UAU^ 1 - ilTVU" 1 , then 

S{Q,A} = S{Q',A'} 

= S{Q + i[R, Q),A} + S{Q, A + i[R, A] - Vi?} - S{Q, A}, (3.9) 

to first order in R. It follows that, for r in the interior of the system, we obtain 

5S 5S 5S 5S 



= -Pm#W ( 3 - 10 ) 

Here we have used a, (3, 7, . . . , for indices running from 1 to 2n (the first n values being 
i, +, i = 1, . . . , n, the remainder i, — , % = 1, . . . , n), with a summation convention, and 
introduced the definition of the current in the NLctM, 



31 



6S\A] 



= (3-11) 



which is implicitly a local function of Q(r) and A(r). The covariant derivative of j CT is defined 
in the same way as that of Q. Therefore, using the notation 



((•■•))= fD[Q):.e?M/Z[A), 



where the • • • represent any functional of Q, and using the fact that the integration measure 
is invariant under U(2ra) gauge transformations, the equation of motion becomes the matrix 
equation, 

D-«r» = 0. (3.12) 

This can be read as the statement that the covariant divergence of the current is zero in 
the interior of the system. As such it is known as a Ward identity, and is a consequence 
of the U(2n) symmetry of the original action S; the existence of a covariantly-conserved 
current as a consequence of a gauge-invariant coupling to a gauge potential is the essential 
content of Noether's theorem. A feature of the NLcxM is that the Ward identities are not just 
consequences of, but are equivalent to, the equations of motion. There are also many similar 
Ward identities when the functional average contains D • j a times other functionals of Q. 
Particular cases of these, including all those that will be of interest in this paper, are those 
that contain other currents j CT . These may be obtained by taking functional derivatives of 
the basic Ward identity ( |3.12| ) with respect to A, since the left-hand side is still a functional 
of A. Functional derivatives of the action yield currents, however, j CT contains A, and so 
does the covariant derivative D, so there are additional ^-function terms. The delta-function 
terms in the response functions that result from the A in j CT will be referred to as contact 
terms, and correspond to those in earlier sections. 

The above Ward identity, or equation of motion, ( |3.12| ), is only the bulk part of the 



system of equations. There are also boundary equations on the reflecting (hard) wall or 
edge. These come from two sources: (i) the boundary terms in the action 5 [A], (ii) the 
boundary term that appears when integrating by parts to transfer the derivative from Vi? 
to SS/ 6 A in eq. ( |3.9| ) when taking the functional derivative with respect to R. The boundary 
part of the equation of motion can be obtained in either of two equivalent ways. One way is 
to take the functional derivative with respect to the full dependence of R on the coordinates 



x and y, and obtain a single equation like ( 3.12 ) but containing delta-function terms at 



the edge. Since an equation that sets a sum of a finite function and a delta function to 
zero implies that each piece separately vanishes, we obtain the bulk equation of motion 
(or Ward identity) as in ( [3.12D , together with a boundary condition that states that the 
function multiplying the delta-function at the edge is zero. The other method is to separate 
the change in the action due to R, after integrating by parts to remove derivatives from 
R, into a bulk part that yields the bulk equation above, and a boundary part, then take a 
functional derivative with respect to the single coordinate x that descibes position on the 
boundary, to yield the boundary conditions. This was the method followed in Sec. IC for 
the quadratic action there, eq. (|1 . 1 1|) . With either method, it is straightforward to obtain 
the results below. 
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First, we should record the actual expression for the current density, obtained from the 
definition ( p.ll ). It is 



#(r) = -iK x Q{x)D»Q{v) - o^DvQdm + io^^Q^y -W)- S(y)}/2 (3.13) 

(we used the identity Qd^Q = —d^Q.Q which follows from Q 2 = I 2n ). Here we see explicitly 
the edge contribution with coefficient cr^' . It is simply proportional to Q at the edge. 

The boundary condition at the reflecting walls (assumed parallel to the x-axis as always 
in this paper), that is obtained by either of the methods described above, can be written 

a° xx QD y Q + a^D x Q + a^°D x Q = 0. (3. 14) 

As in the case of the bulk equation of motion, in the quantum field theory of the NLcrM 
(defined by the functional integral), this is valid only when inserted in the average ((. . .)). 
This equation can be interpreted as stating the covariant conservation of the current at the 
edge, in a very similar way to that discussed in Appendix A within SCBA. For cr X y'° = 0, we 
would have only the first two terms, and it would state simply that the normal component 
of the (bulk) current in eq. ( |3.13D tends to zero at the edge. These terms originate from the 
edge term left after integrating the bulk part of the action by parts in order to take 5 /SR. 
In the presence of a non-zero &iy , this boundary condition is modified to include the last 
term, which is the covariant derivative (along the edge) of the edge part of the current in 
eq. ( |3.13|) . Thus eq. ( |3.14|) is equivalent to 



.7 y, bulk Acj^edge — 0> (3.15) 

where the edge contribution is obtained as j Xte ^ se (x, W) = Jw-o+ dy j x (x, y), and similarly 
for the edge at y = 0, as in Appendix A. The edge term originates, of course, from taking 
5/5R on the edge term in the action itself, eq. (|3.8| ). In the boundary condition, the last 
two terms can be combined to leave 

a° xx QD y Q + a° xy D x Q = 0. (3.16) 

This is the analogue of the tilted boundary condition discussed in Sees. I and II, generalized 
to the full nonlinear field theory, and including A; as in those discussions, only the total 
a xy enters the boundary condition. We will show, in Sec. IIIC below, that, in leading order 
in perturbation theory, this boundary condition reduces to exactly the one used in earlier 
sections. 

To close this subsection, we obtain Ward identities that apply to moments of the two- 
probe conductance, in which the currents are integrated across sections parallel to the y 
axis. We have already seen that these should be calculable using only Green's functions at 
the Fermi energy, which can be obtained using our generating function or the NLcrM. As in 
Sec. IIIA, we will therefore here specialise A to be in the x direction and independent of 
y. Functional derivatives of Z[A] of the form S/SA x (x) then produce the desired mean and 
variance of the two-probe conductance. If we specialise to such A in the NLcrM action S[A] 
(noting that A x (x, W) = A x (x,0) = A x (x,y) for all x, y), then we see that the a xy terms 
can be combined, using an integration by parts, to leave 
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S[A X ] = J d 2 r [~<T° xx ti {D^QD^Q) - ±a° xy tr {t, v QD ,QD V Q) + V ti (QA) 



(3.17) 



in which, we emphasize, A takes the specialised form A(r) = (A x (x),0), independent of y. 
This action is gauge invariant, and A x remains y-independent, for gauge transformations 
U that are independent of y. Using such a transformation, we can obtain, similarly to the 
above derivation, the identity 

D x Jdy((f x (x,y))) = 0, (3.18) 

in which j x and the action SLA^] to be used in calculating the average still contain A x . We 
note that, in J dyj x , using (|3.13| ), cr x 'y and cr^'° terms can be combined into a single term, 
as implied by the action S[Ae], in which the same is true. Thus from this point on, only the 
total (7° , enters the calculations for the conductance and its moments. 

y 



On taking further functional derivatives 5/5A x (x) of ( 3.18 ), we can obtain identities 
involving the mean and variance (or alternatively, second moment) of the conductance. 
From the discussion in earlier sections, we should have 



^(1,^ = 0, (3.19) 

for all x, x 1 inside the sample, including x = x', which expresses the independence of the 
conductance on the location of the cross-sections, and similar statements should hold for the 
higher moments and for the dependence on the other variables x' , .... One may be concerned 
that D x appears in ( |3.18|) , not d x , and that this might lead to additional terms containing 
e.g. 5(x — x') when further functional derivatives are taken. However, we have verified 
that such terms vanish, to all orders in perturbation theory, for the replica components we 
require to produce the mean and variance of the conductance, as in (|3.4j ) and Q3.5|) . Thus, 
we have performed all the steps in a derivation showing that the calculation of the two-probe 
conductance and its variance, in the rectangular geometry, can be carried out within the 
NLaM, including the independence of the locations of the cross-sections, which allows the 
use of an average over these locations. The explicit expressions for the mean and variance 
are given in Sees. IV and V, and evaluated to leading order in the perturbation expansion. 



C. Perturbation expansion, current conservation, and bilocal conductivity 

The approach given in the previous subsection is sufficient for the two-probe conductance, 
provided the cross-sections used are parallel to the y axis. For more general cross-sections, 
and to consider the non-Fermi-energy effects and the bilocal conductivity, a deeper analysis 
is required, which is contained in the present subsection, but which may be skipped (apart 
from the first part introducing the perturbation expansion) by readers interested only in the 
two-probe conductance. 

Our goal here is to compare the consequences of the definition of the current and the 
Ward identities with the properties of the currents in our model of the original electron 
system. To provide motivation, we will compare results in the NLaM formulation with 
those in the SCBA in the previous section, and for this purpose we will now introduce the 
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perturbation expansion of the model. We then show how to modify the coupling of A to 
the NLcrM so as to reproduce the properties found in Sec. II. 

The matrix Q, which obeys = Q, Q 2 = l2„, can be parametrized in the following way: 



M * -vM- (320) 

where z is an n x n complex matrix. Expanding the NLcrM action in terms of the ^-matrix, 
we get 

S[A] = S [A] + Si [A]. (3.21) 

where So [A] is the part quadratic in z, z^ and A, which is the same as in eq. ( |1.11| ), except 

that there are now n 2 copies of z, and we include the gauge potential A. For A of the 
restricted form 

a =(a°+ A o")- < 3 - 22 > 

where A +_ is a complex n x n-matrix-valued vector field, and A h is its adjoint (these are 
the only components of A that will be used below), we have 

rO 



S [A] = - a ~fJ d 2 rtv ((d,z - 2iA+~)(d^ + 2iA~ 

j n 

J dVe^tr ({d^z - 2iA+-)(d u z jf + 2iA; 

n • 11,0 

J dVe^tr^cU) - % ^-jdl lx tT{A+-z'< - zA~+). (3.23) 
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The line integral § dl^ is taken in the counterclockwise direction around the edge of the 
sample. Here and below, we use the symbol tr for a trace on the n-dimensional space, 
as well as for that on the 2n-dimensional one; it should be clear from the context which 
is meant. Si describes the interaction between the diffusion modes caused by quantum 
interference effects. We give it here only for A = 0, and to the order 0[(zz^) 3 ] required for 
our later calculations, 

Si[0] = I d 2 r ^-^ tl [d^)d^) + d,(zh)d^z)] 
^-tile^d^zz* zd v z ] ) - e^d^(z ] 'zz x d v z)\ 
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ti[d^(zz ] )d^zz ] zz ] ) + d^(z ] ' z)d^(z ] ' zz ] ' z)] 



CT 

" xx 



64 

-^tvie^zzhzhd^) - €^d^zzhz^d„z)} + 0[(^) 4 ]| . (3.24) 

Notice that terms proportional to erJL can all be written as total derivatives, therefore they 
can be expressed as boundary terms. To calculate the ensemble average of any quantity 
X[z, z*], we perform the following expansion: 
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r 00 1 

((X)) = lim / D[z,zt)I[z,z*)X[z,zi]e*>W» £ -p^M]. (3.25) 

Here I[z, z*] is the Jacobian needed to make the measure in the z, z* space invariant under a 
U(2n) rotation at each r. The explicit form of this Jacobian will not be needed. Its only role 
is to cancel quadratically-divergent diagrams that arise in perturbation theory, in a manner 
that is standard for all NLaM's (see, e.g., Ref. [0]). The terms in the expansion can be 
written in terms of averages calculated using the quadratic action with A = 0, defined by 

«•■■»<, = jD[z,zt].--e*>^/Z [0], (3.26) 

in which the functional integral in the denominator is the same as the numerator but with 
the insertion • ■ ■ omitted. The expressions can be evaluated by contracting pairs of z and 
z\ which gives the diffusion propagator, 



^«4(r)z„(0»o = <W(r,r'), (3.27) 



which obeys the same conditions ( |1.16| ) as in earlier sections. The basic perturbation ex- 
pansion is now a series in powers of l/cr° x , though it will also be convenient to expand in 
powers of 7 = cr° /cr^., to obtain a double expansion. 

We now return to the physical meaning of the Ward identities (|3.12j ), ( ft.l6|) , that resulted 



from the gauge invariance of the action S[A\. We wish to compare these with our physical 
expectation that the current is divergenceless inside the sample, and that no current flows 
in or out of the sample at the reflecting walls (as we have shown in Sec. II, and discussed 
in Sec. I, the current response obtained in the SCBA is not divergenceless, because of the 
bulk &xy'° term, but that is a non-Fermi-energy effect that will not be considered in the 
present formalism until later in this subsection). The first difficulty that seems to arise (as 
with eq. ( |3.18|) ) is that the Ward identity states D • j CT = 0, not V ■ j = 0, as we might have 



expected. The vector potential A present in D will generate 5-function terms when further 
functional derivatives are taken to obtain Ward identities, as must be done for the bilocal 
conductivity tensor and analogous correlators of more than two currents. However, as we 
mentioned in connection with the conductance in the previous subsection, in practise, for 
the particular components of A that yield the physically-relevant conductivities, this does 
not seem to occur. For example, in addition to the results cited in the previous subsection, 
we can show that the Ward identity implies 

d„{a^(r,r')) = (3.28) 

for the mean bilocal conductivity tensor calculated in the NLaM using the action S*[A], and 
this is valid for all r, r' inside the sample, including r = r', to all orders in perturbation 
theory. 

The boundary condition (|3.14j) also involves the tangential covariant derivative of the 
edge current, not the usual partial derivative as one would want in the electronic system. In 
this case, we do find a clash between the theory as formulated and our intuition. Eq. (|3.15|) 
is more explicitly 
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jy,bulk ®xjx,edge ~ O ^ X ' ^ ' (3.29) 

The left-hand side is the combination one might have expected to be zero. However, it is 
nonzero when A x is nonzero, implying that conservation is violated by 5-function terms on 
the edge in the bilocal conductivity and its moments. Note that similar commutators [A x , Q] 
and [A y , Q] appear in j^ bulk but are not a problem. 

To illuminate the point further, we can use the perturbation expansion and compute the 
mean bilocal conductivity within the NLcrM as formulated so far. From So [A], the equation 
of motion for z' (formulas for z are similar) is, in the bulk 

alVV = -2i (a° xx d,A- + + , (3.30) 

(as in Sec. IA, where, however, a xy = a x ^ and a X y'° = 0), and at the edge is the tilted 
boundary condition 

<j° xx (d y J + 2iA-+) - al y {dj + 2iA~ + ) = 0. (3.31) 

These equations give the generalization of the "classical" theory of Sec. IA to include the 
edge currents with coefficient cr^'°, on using the identifications given in Sec. ID. The currents 
are given, in the present approximation, by 

6S /6A*T = j~ + = ~z[<4(^+ + 2* V) + <4V(^ + 2iA' + ) 

+ a I x I y °(5(y-W)-S(y))S^]. (3.32) 

The bulk equation of motion can therefore be written d^j~ + = 0, while at the edge we have 

iy,bulk — ®xjx,tdge = ~~ a xy A x + , (3.33) 

which can also be written 

Jj/,bulk ~~ ^ix.edge = 0' (3.34) 

where the covariant derivative is the linearized version of that in the full NLaM, namely 
= + 2iA~ + , while D fl D u z J< = d^D v z^ . Thus, for current j~ + as defined here, 
current conservation in the naive form is violated by 5-function terms at the edge. This can 
be rectified, but before doing so, we calculate the bilocal conductivity of the present model 
in the same approximation. 

The mean bilocal conductivity in the present approximation is obtained as (compare eq. 
( |3.4| ); we leave implicit the choice of all replica components equal to 1, and the n — > limit) 

(a^{r,r')) = - Urn M _ +(r/) (0,7 + W))o 

= Kx^ v + k' y \Mr ~ r') - (0; + (r)j+-(r'))) . (3.35) 

On evaluating this using eq. ( |3.32| ) (with A = 0), the similar formula for j + ~ , and the 
definition ( |3.27|) of d, we obtain the same result as in eq. ( |2.20| ), except that the bulk <t^'° 
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term is not present. This result obeys d^a^r, r'))o = for all r, r' in the bulk, and eq. 
( 13.33| ) implies that 

rW+0+ 

(cr yi/ (r, r')) -d x dy (a xv (r, r')) = a I x Ifi 5 ux S(r' - r) (3.36) 

Jw-o+ y 

for r at the upper edge y = W (and similarly for the lower). In effect, in the edge channel, 
A~ + simply creates current, so the naive conservation law is violated. 

In the SCBA, we did not directly address this issue, but derived d x j x ^^, ~ Jy,buik = f° r 
r^r' only. On the other hand, in the SCBA we also found a non-Fermi-energy contribution 
<j I x I y )0 e^ u 5{r — r') in the bulk, which means that in the presence of E, there is an additional 
part [e 1 jb^a^e^Ey in the bulk current 

■SCBA _ -SCBA , e _//,0. ip (n o 7 \ 

Jfj, — Jfi,E=E F ~r %y e ^" C/ ^' yo.ot ) 

within the SCBA. If we introduce a corresponding change in the bulk current here, so that 

J^,mod = + ~ a xy° e ^^v + 1 (3.38) 

then the tilted boundary condition Q3.33j ) becomes 

i?/,bulk,mod — ^ia;,edge,mod = (3.39) 

in the presence of A. Thus the modified current is conserved (not covariantly) at the edge 
(and so is j^ CBA ); the current in the edge channel comes from the bulk. Physically, there 
are two modes of conduction response to an electric field in the system. One is the "sliding" 
of the total charge density, which gives the bulk Hall conductivity cr^'°- This is a non- 
Fermi-energy effect, and is a local (5-function) response to an electric field. The other is 
the Fermi-energy response, which is diffusive in the bulk (including the Hall effect with 
coefficient a x 'y) and is chiral along the edge. As discussed in Sec. I, the a X y'° bulk effect 
implies V ■ j 7^ 0, meaning that dp/dt 7^ 0. At the edge, there is no charge accumulated. A 
tangential electric field at the edge can produce a bulk current normal to the edge, and also 
a Fermi-energy edge current that increases along the edge. These effects involve the same 
coefficent <7 X y'°, and the result is that no current is created, so no charge accumulates at the 
edge. This occurs because of a version of the Laughlin-Halperin gauge-invariance argument 
pILEDI] . A change in the potential (which is essentially what z is) would accumulate a charge 
density of order the inverse velocity of the edge states, but the same velocity also appears 
in the edge current, which carries away the charge. 

We now propose a modification of the NLcxM action which incorporates this non-Fermi- 
energy effect so as to recover the SCBA bilocal conductivity tensor in full, and maintain 
current conservation at the edge (though not in the bulk), to all orders in perturbation 
theory. Our proposed action (in which we reinstate the full In x In matrix A) is 

5 mod [A] = J d 2 r {~olMD»QD,Q) - l -a^tv{e^QD^QD u Q) 

~ IvilM^vQd.QdyQ) + ^J°tr(e M „g[i M , Q] [A u , Q}) + 7/tr(QA) 

+ / dx {tr(A x (x, W)Q(x, W)) - tr(A x (x, 0)Q(x, 0))} . (3.40) 
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The added a^y term maintains SU(2n) global, but not local gauge, symmetry, corresponding 
to the nonconservation of the corresponding modified current which contains an additional 
term in the bulk: 

J/I,mod = = J ft ~ 2 aI Jy° e lJ,v\Avi Q\- (3-41) 

The modified equation of motion must be obtained from SS mo< i{Q', A}/5R = (as in Sec. 
IIIB) without using gauge invariance. It can be written in terms of the modified current, to 
give the modified Ward identity in the bulk, 

A^mod = V R QA,, Q}) + i[A„ [A v , Q))} . (3.42) 

The modified current is not covariantly conserved, because the modified action is not gauge 
invariant. However, the boundary condition eq. ( |3.16D is unchanged, because the added term 



in S mo d contains no derivatives, so does not give rise to any boundary terms. Nonetheless, 
the interpretation of the boundary condition changes, because the current has been modified. 
In terms of the modified current, the boundary condition states that at the edge 

jy,mod ~ ^xi x , edge, mod — 0' (3.43) 

which is "current conservation" . Strictly, our arguments imply that this modification of the 
action applies only for the H — and — h components of A; for the other components, the 
correct form may depend on what is assumed in the underlying model the NLctM is supposed 
to represent. 

For the quadratic part So,mod of the modified action S mo d, the corresponding formulas 
are: for the action, 



S , m od = So - << u J ^re M ,tr(A+-A; + ); (3.44) 
for the current, 

i/^,mod = Jfi ~ °xy ^^v^ u , (3.45) 

as in eq. ( |3.38| ) above; for the modified equation of motion, 

d^mod = -< (3-46) 

which corresponds to the earlier eq. (|1.27|) ; and in the bilocal conductivity tensor, the bulk 
contact term cr^'°e Atl/ 5(r — r') appears as in the full SCBA result, eq. fl2.20|) . The boundary 
conditions on z, are, however, unmodified. 

When using the modified action for calculations of conductance, there is no change to 
the results, as long as one uses cross-sections that are parallel to the y-axis, and therefore 
the expressions contain the integrals of the ^-components of the currents, as we were doing 
earlier. For those calculations, we already showed at the end of Sec. IIIB that the conduc- 
tance is independent of the positions of the sections. Since the vector potentials used there 
have x-components only, the extra term in S mo( i is zero. In addition to the usefulness of the 
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general S mo d for reproducing the bilocal conductivity tensor and maintaining the current- 
conservation properties at the edge to all orders in perturbation theory, it is also crucial 
for the conductance if one calculates the flux of current through more general cross-sections 
than those specified above. In general, a cross-section could be any curve that intersects 
the edges of the sample just twice, once on each of the reflecting walls. Two such sections 
may intersect at isolated points (instead of along their whole length), and the intersections 
are then said to be transversal, that is, the normals to the curves at the point of intersec- 
tion are not parallel (nor antiparallel). The question arises whether the conductance and 
its moments calculated using such sections is independent of the position and shape of the 
sections, including the case in which they intersect. When they intersect transversally, the 
bulk crly contact term will contribute at the intersection point, even within the SCBA. 
Our preliminary investigations of this, which will not be included here, show that these 
contributions are needed to cancel effects of the a™' terms at the edge, so as to maintain 
conservation of the total current, and that the conductance obtained is the same as for the 
straight sections, for any shape and position. We believe this to be true in general, to all 
orders in perturbation theory. This shows that the use of the modified action is obligatory 
for such general calculations. Note that, for more general geometries, such as four probes, 
such intersecting sections will be common. 



D. Further details of perturbative calculations 



1. Boundary perturbation expansion for the diffusion propagator 



It is difficult to calculate the diffusion propagator as defined by eq. ( P-.16Q explicitly 
(however, see reference fl39|), although propagators for simpler geometry such as a half- 
plane f4"B|| , an infinite strip ||8],|33] or an annulus |58| can be found (the results for the 
infinite strip can be obtained by conformal mapping from the half-plane |48|). We can 
perform a boundary perturbation expansion in powers of 7 using the propagator at 
7 = 0, which we define as d°(r, r'). d°(r, r') can be constructed out of the solutions for the 
following eigenvalue problem: 



- V 2 0° = A0° 



(3.47) 



with the boundary conditions d y (f) (x,W) = d y <f) (x,0) = and <p°(0,y) = (f)°(L,y) 



0. The eigenfunctions are 4>^ m (x,y) 



eigenvalues A° r 



ILW 



cos(niry/W)sm(m'Kx/L) with corresponding 



(wr/VF) 2 + (rrnr/L) 2 , where n 

4>nm( x ,y) ( t>'nm( x ' \v') 



0, 1, 2, m = 1, 2, .... We have 



Ar,r') = EE 

n=0 m=l 



A 



(3.48) 



Using the bulk diffusion equations for d°(r, r') and d(r, r'), we get 



d(r 2 ,n) =rf°( ri ,r 2 )+ / dS- 

Jc 



d°{r, r 2 ) W(r, n) - d{r, n) W°(r, r 2/ 



(3.49) 
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where C is a closed surface enclosing the disordered region including the edges (see figure 
[l]). Let us divide the surface C into four parts C%, C 2 , C 3 and C 4 . Applying the boundary 
conditions for the propagators for different sections, we have the following: 

d(r 2 ,r 1 ) = d°(r 1 ,r 2 )+ / dS n d°(r, r 2 )9 n d(r, n). 

Plugging in the boundary condition d n d = jd^d, we have 

d{r 2 , n) = d°(r a , r 2 ) + 7 / d°(r, r 2 )<9 t d(r, n). (3.50) 

The above equation generates an expansion in terms of d° and in powers of 7. Using B to 
denote the operation Jci+c 3 dS n jdt, we can write schematically: 

d = d° + d°Bd° + d°Bd°Bd° + .... (3.51) 



2. The ID propagator 

The above expansion is not valid for extremely narrow samples with W <C L, where 
equation ( |3.51| ) needs to be summed exactly. However, the propagator itself approaches the 
ID diffusion propagator with conductivity (J XX {1 + ^ 2 )W . Such limiting behavior can be 
demonstrated by rewriting the quadratic action in the following way: 

*(,, = *r {(1 + fM M + trKV + *W( V - 0*,')]} . (3.52) 

In the limit W ^ L, the main contribution to ((z Jf z)) comes from the low-lying eigenmodes 
of —V 2 which satisfy (d y + ■yd x )4> L (r) = 0, (d y — •yd x )4> R (r) = in the entire strip. The 
eigenvalues of these modes are separated from those of the other modes by a gap of order 
L 2 /W 2 P5| . In other words, the term / d 2 r tr[(d y z + r yd x z){d y z^ — •yd x z ji )} in the action can 
be ignored in the ID limit: 



Jim J D[z,J]zy kl e s ° = J ^D[.,,t]4.^ e -¥(l+7 2 )W/dxt r [^.t] 

= 45 a 8 jk d 1D /a° xx W : (3.53) 

where 

= 4 ~ sin(m^)sin(^x7L) ^ 
(l+7 2 )L^r 1 (mix /L) 

More generally, by symmetry, we obtain the NLaM in one dimension (where no topological 
term is possible) with cr x ' x D = cr xx (l + ^ 2 )W = W/p xx as the coefficient in the action. 
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IV. TWO-PROBE CONDUCTANCE 



A. The boundary contribution 



From the nonlinear a-model, the average conductance is 

5S 



(d) = -75 . !i n m n [ dx 1 [ dx\ ({ jj+zry 

L z A-^0,n-»OJo JO oA xll {X\ 



5 2 S 



7T», (4- 1 ) 



where 



6S 



5 2 S 



A=0 



w 



dyi 



dAt^x^SA-^x^) 



fW 1 1 1 

-J dyxalJixx - -{zz ] )u - -{z ] z)u - — (zzW) 



11 



~(zWz)n + iWuWu + 0(z 6 )]. 
10 o 



(4.2) 



To leading order in l/cr^, we get 

S ° = 1? - r') - +7W, - 7^)«4i(r)^n(r')))o] 



d 2r ' °lx 



5(r - r') - (d x + ld y ){& x - 7 ^)d(r, r') 



(4.3) 



We have thus recovered the result from the diagrammatic expansion, and the other ap- 
proaches described in Sec. I. Therefore the following results apply to any of these approaches. 

The long-ranged term in the above expression comes from the ladder diagrams (i.e. 
diffusion) and in the absence of magnetic field, it does not contribute to the two-probe 
conductance . In the presence of magnetic field, this is no longer the case. One can show 
that the local term gives the Ohmic conductance a xx W/L. The long-ranged term, which 
involves volume integrals of total derivatives, gives the difference of boundary values of the 
diffusion propagator. Since the diffusion propagator goes to zero in the leads, the d x and & x 
terms vanish upon volume integral. We are left with the boundary difference at the upper 
and lower edges, 



y xx 1 w xx 1 



L pL 

dx / dx' [d(x, W\ x', W) + d(x, 0; x', 0) 
Jo 

-d(x, 0; x', W) - d(x, W; xf, 0)] . 



(4.4) 



The (7° -dependent part is expressed as a boundary term and vanishes when the magnetic 



field is zero (7 = 0), or if the system is subject to periodic boundary condition in the 
transverse direction. 

At one-loop level, which is the next order in l/cr xx , we have verified that the interference 
correction to (g) vanishes in the limit n — > 0. Therefore the presence of edges does not 
change the conclusion of the previous perturbative calculations P^-j27| that there is no 
weak localization correction to a xx of relative order l/cr xx . In general, we do not expect the 
presence of edges to have any effect on the renormalization flow of a xx in the perturbative 
regime, since it is dominated by short distance effects in the bulk. Whether a xy is ever 
renormalized perturbatively when the system has edges is less clear to us. 
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B. The small 7 correction for the two-probe conductance 

For small 7, we can make use of the propagator d° at 7 = to obtain the leading 
correction to g°. Plugging d°(r, r') in the boundary term of the two-probe conductance, we 
get 



9\l) = a°„W/L[l + l 2 h{L/W)} + 0( 7 4 ), (4.5) 



where 



Ri ,2 oo 00 , , f U((3)L/k"W, W > L; 

EE 33-10- 1/2, W = L 



71-4 ^ 2 m=l,oddn=l,odd m2 m2 + ^2"- | ]_ W <^ L. 

(Here £(s) = Em=i m ~ s is the Riemann zeta function, and we note that £(3) ~ 1.202.) 
Thus, for an extremely wide sample, W 3> L, the effect of the edges can be ignored. In the 
ID limit (W/L — > 0), we get g° — > (1 + 7 2 )<t° ;e W^/L, which is consistent with our result that 
in the NLaM, <7^ X D = (1 + , ~f 2 )cr xx W. This is also consistent with other results, as noted in 
Sec. IA, valid for arbitrary 7, which show that the mean conductance in the ID limit can 
be obtained from this one-dimensional conductivity cr®£ D = W/p xx . 



V. VARIANCE OF THE CONDUCTANCE 



In this section, we evaluate the variance of conductance to leading order in 1 
the nonlinear a-model, we get 

5 2 S 5 2 S 

( ( 5A+- n ( Xl ) (4 ) 5At; 22 (x 2 )5A-+ 2 (x' 2 ) 
5 2 S 5 2 S 



/a° xx . From 



where 



u o o 

5At M (x! ) 5A~% 2 (4 ) SA+ 22 (x 2 ) 5A X ^ (x[ 
5 2 S 5 2 S 



5At; u { Xl )5A + x ; 22 {x 2 ) 5A X +, WSA-^ix't) 
5 2 S 5S 5S 



' SA+nixJSA'+M) 5At; 22 {x 2 ) 5A x + 2 (x' 2 ) + 1 ^ 2 

S 2 S 6S 5S 

' SAi^xJSA-^x',) 5A&(x 2 ) SA-^x',) + 



5 2 S 



SS SS 

'5A+n(xi)5A+2 2 (x 2 ) ^,22(4) 5A^(A) +{+ ^~ } 
SS SS SS SS 



OD OD OD OD 

" m5R sa^m) 5At; 22 {x 2 ) 5A-; 22 (x' 2 ) } } — tcd 



(5.1) 
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5 2 S r w a 

L ^^-^^44, (5.2) 

and the other functional derivatives are all evaluated at A = 0. The leading diagrams are 
shown in Figure f|. Various vertices are denoted by polygons with the wavy tail indicating 
while the lines linking the vertices are the diffusion propagators. The diamond shaped 
vertex with no wavy tails comes from the 4-point interaction term in Si. We have also 
obtained the same set of diagrams using the diagrammatic approach. The diagrammatic 
approach is complicated because various vertices need to be evaluated separatedly. In the 
presence of magnetic field, the vertices are dressed with non- vanishing G + G + and G~G~ 
ladders, although in the end they can all be expressed in terms of a xx and a xy . For the 
NLaM, the vertices can be obtained from the action. Figure [5] shows that one particular 

vertex from <- A+ - c *-+ is equal to the sum of four diagrams in the diagrammatic approach. 
Figs, ^a and |4]b are the only diagrams considered in previous UCF theories. The rest of 



the diagrams have been considered by Kane et al. |38[ and it is known that, for 7 = 0, they 
give rise to the long-ranged correlation in local current response but they do not contribute to 
the variance of conductance in the absence of magnetic field [ . One can show that these 



additional diagrams can all be written as boundary contributions and they vanish when 
7 = for the same reason as the ladder series vanish in the case of the average conductance, 
when written in the area-averaged form. However, in the presence of the magnetic field the 
additional diagrams give rise to Hall-ratio-dependent contributions. The work of KY and 
ML discussed the effect of the tilted boundary condition on the diffusion propagator but did 
not consider the additional diagrams. 
From Figs. f|a and |]b alone, we get 

(Sg 2 ) a ,b = y± {4Tr(^ T ) + 2Tr{dd)} . (5.3) 

Using the classical network model of Sec. IB, we have calculated the diffusion propagator 
d for a range of values of 7 and W/L. Using this propagator, we find that Tr(dd) and 
Tr(dd T ) are smooth functions of 7 and L/W. The peaks reported in ML in the variance of 
the conductance (as given by eq. (|5.3p ) are not observed in our exact numerical calculation. 
The argument advanced by ML for the existence of "resonant" peaks due to the tilted 
boundary condition effects is not supported by this calculation. We emphasize again that, 
in any case, eq. Q5.3p is not the full expression for the variance of the conductance, because 
there are other diagrams that were omitted by KY and ML. 



A. The recovery of UCF result in the ID limit 

The importance of the additional diagrams can be best demonstrated in the quasi- ID 
limit (W <C L), where a xx and a xy combine to form a single parameter, (J XX D - This limit 
is well described by the random matrix theory of the unitary ensemble. For general reasons 
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given earlier, we expect the variance of the conductance to approach the well-known ID 
UCF result, independent of the value of 7. 

Plugging in the ID diffusion propagator of eq. ( |3.54|) , we get 



where we used 

This is essentially the argument used by KY and ML, except that they gave versions appli- 
cable at finite temperature. However, this result of these authors, that the variance of the 
conductance depends on 7 even in the ID limit, is incorrect. 

Among the additional diagrams, Figs. |]c, d and the sum of |]i and i' vanish to order 7 2 
for all W/L; Figs. |]g, h, i, i', j and j' vanish as W/L — > 0. The dominant contributions come 
from Figs. |]e and f: 



m«f*D 27^/1(00) = -( 2 y). (5.5) 

Figs. §e and f thus cancel the j 2 correction from Fig. §a and b. We get in total, to order 7 2 , 

($g 2 )w/L^o = jg- (5.6) 

Thus, in ID the UCF result of 0,0,0] is recovered, at least to order 7 2 . We remind the 
reader that the ID UCF result holds only when the length L is less than the ID localization 
length, £i D , and that £ 1D is of order a^ x D = Wj p xx , which is much larger than the lower 
limit (W) on L in the diffusive regime p° xx <^ 1. 



B. The variance of the conductance in 2D 

For a wide sample with W/L arbitrary, the variance of the conductance depends on 
the Hall ratio. We will calculate the correction to the usual result for the 7 = unitary 
ensemble, to order 7 2 . In the 2D limit, the individual diagrams, Figs. |]e-h, j and j', all have 
logarithmically-divergent parts, however, their logarithmic contributions cancel out. (The 
cancellation is guaranteed by the fact that S is not renormalized at one- loop level.) There 
can be even more divergent diagrams containing 5(0), but these diagrams are canceled by 
diagrams generated by the measure I[z, z*]. (Since they are at least of order 7 4 , they are 
not explicitly calculated in this article.) The total 7 2 correction is finite. Summing up the 
contributions from Figs. p§a-j', we get for a square sample 



I _ 2 256 



9.06— + 2A0Y 



IT 4 ' 7T 8 



+ 0( 7 4 ). (5.7) 



The expression for the variance for arbitrary W/L is given in Appendix D. 
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VI. CONCLUSION 



In this paper, we considered the mesoscopic conductance and its fluctuations in the 
presence of a magnetic field for a realistic two-probe geometry. Our perturbation theory has a 
different structure from previous theories P,|7|,|32| because of the presence of two conductivity 
parameters, a xx and a xy . We found that a xy not only enters the boundary condition for 
diffusion, as was noted in Refs. f33H35l, but also appears in the current vertex and other 
vertices which govern the interference processes. As a result the two-probe conductance 
and its variance in the perturbative regime depend on the Hall ratio 7 = c xy /a xx . ® m 
calculations differ from the previous results P5| , |54]| since we have not only modified the 
boundary condition but also considered additional diagrams which vanish in the zero field 
limit or in an edgeless system. Our main result is that the UCF are modified in the presence 
of edges; the variance of the two-probe conductance, although it is still of order 1, increases 
with the Hall ratio, as shown in eq. ( |5.7| ). However, in the quasi-lD limit of a long sample, 
the usual universal result is recovered. 

The reflecting boundary condition at the "hard" wall (or "edge" ) is crucial for the depen- 
dence of the conductance on the Hall ratio 7 that we find. If this is replaced by a periodic 
transverse boundary condition (i.e. a system on the surface of a cylinder), the results of the 
usual unitary ensemble in 2D are obtained; the results of Xiong and Stone are easily 
modified for this case, for which they are correct. While a cylinder may seem hard to realize 
experimentally, it can be mapped to an annulus by a conformal mapping. The annulus is 
sometimes known as the Corbino disk, in which there are no edges, and a radial voltage drop 
is applied to induce a current flow. Thus for the disk, the conductance fluctuations should 
be a universal function of the ratio of the inner and outer radii, with no dependence on 7. 

The experimental observation of the effects we find depends first on being in the regime 
L in 3> L, W, so the system is phase coherent, and on having an elastic mean free path I 
due to impurities such that L, W ^> I. Our calculations only address the metallic regime 
of conductance fluctuations at large diagonal conductivity a xx , where perturbation theory 
is valid. In principle, this approach is valid for any value of the Hall ratio 7 = o~ xy /a xx , or 
of the Hall angle 9h = tan -1 7. For simplicity, we expanded most of our results also to first 
nontrivial order in r j 2 . The terms in (l/cr^,) 2 that are left out cannot be neglected if the 
system size L or W exceeds the order of £, pe rt, the crossover scale at which the renormalized 
conductivity becomes of order 1 or less. If L, W are greater than £, per t, the system crosses 
over either to the localized regime where a xy becomes quantized, or, for Fermi energies near 
the critical values that lie near the centers of the Landau bands, to the critical transition 
region between the plateaus; our theory does not apply to either of these. Therefore, one 
must use mesoscopic systems that are not too large. Fortunately, since £ pert ~ le^^ 2 , this 
is not difficult if a xx 3> 1. According to the SCBA results reviewed in Sec. IIB, a xx will be 
large unless either the Landau level index iV of the highest partially-occupied Landau level 
is of order 1, or the Fermi energy lies in the tail of the density of states of the disorder- 
broadened Landau bands, when u c r Q is large enough that these are well developed. Thus the 
magnetic field B must be large enough to suppress the Cooperons, so the system is in the 
unitary (broken time-reversal symmetry) regime, but not too large. (We do not generally 
require u c t > 1, though this would ensure that 7 > 0(1).) In effect, for the observation 
of the effects found in our theory, ideal conditions would be that the system should exhibit 
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Shubnikhov-de Haas (SdH) oscillations, but not well-developed quantized Hall plateaus, even 
for asymptotically low temperatures. As the Fermi energy or magnetic field varies through 
a Landau band, yielding such an oscillation in a xx , a xy varies monotonically, which implies 
that the ratio 7 = cr xy /a xx varies. There is therefore a lot of scope for varying 7 by varying 
either the field B from low values (7 ~ 0) to larger, or as the field sweeps through a single 
SdH oscillation. However, since the amplitude of the fluctuations depend on 7, it will be 
necessary to collect statistically-independent values of the conductance without changing 
7 too much. Thus the simplest experimental method, which uses magnetic field as the 
ergodic parameter, will not work and some other technique must be used to vary the sample 
conductance at fixed B. Finally, as the quantized Hall plateas are reached, localization 
effects will suppress fluctuations strongly between the centers of the LLs, and our theory is 
not applicable (although such measurements would be interesting). 

While the calculations in this paper have addressed only the weak-coupling regime at 
a xx 3> 1, it is interesting to speculate about the effects of a xy on the conductance and 
its fluctuations in the critical regime of the integer quantum Hall effect, when the system 
has edges. The critical regime can be defined by the conditions Lj n 3> L, W, and L and 
W between £ pert and £, where £ > £ peri is the localization length, which diverges as Ep 
approaches any of the critical values E cN , N = 0, 1, . . . . We expect that the renormalized 
local conductivity parameters a xx , a xy , a xy are still meaningful, and that a xx and a xy = 
a xy + a xy t & k e on universal values (=1/2 (mod 1), in the case of a xy ) at the critical points. 
This raises the question of the renormalization of the two pieces a xy and a xy , and whether 
the values of these are universal separately at the critical point. We note that at the localized 
fixed point, the behavior may be described by saying a xy = 0, a xy = (mod 1), so that these 
parameters do approach universal values in this regime. For a xx , there is a widespread belief 
that it takes the universal value 1/2 at the critical fixed point, though it is not always clear 
if the calculations done to support this are describing the local conductivity parameter a xx , 
rather than a mean conductance in a particular geometry. The relation of these is not known 
in the critical regime at the present time, and, as we have seen, is not simple even in the 
perturbative regime, if the system has edges. We expect that for a two-probe system with a 
periodic transverse boundary condition, a xy should not contribute to the conductance in the 
critical regime, just as it does not in the perturbative theory in this paper. Even then, the 
mean conductance is not in general given by a xx W/L, since non-Ohmic behavior is expected 
at least for L>If where the system approaches a localized quasi-lD limit. Thus, even in 
the case of a square sample with W = L and periodic transverse boundary condition, it is not 
clear that (g) = a xx . The effect of the edges in the critical region is nicely shown in a recent 
paper [|23[| , which examined the mean, variance, and distribution of the conductance in a 
two-probe geometry like ours, with W = L, and for both reflecting and periodic transverse 
boundary conditions, i.e. with and without edges. The results show that the boundary 
conditions do make a difference (however, finite size effects are significant, as shown for the 
periodic transverse boundary condition case in Ref. [p2]|). The authors tentatively attribute 
this to "edge currents," but as we have seen in the perturbative regime, there are edge effects 
(described by a xy ), that are not solely due to edge currents carried by edge states (which 
are described by cr xy ). The boundary effects make themselves felt throughout the system, 
due to the long-range correlations in the critical regime. They are relatively unimportant 
only when W ^> L. In fact, dependence on the boundary conditions, say on whether they 
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are periodic or reflecting, would occur even in the absence of a xy , as it does in the weak 
coupling regime (see e.g. f32|). A further implication, suggested by our results, is that the 



critical conductance properties, in a given geometry that possesses edges, may depend on 
the fixed point value of a xy , i.e. on which transition is being studied. While the structure 
of the critical field theory [including a xy (mod 1), a xx , and the critical exponents] should be 
universal, this may not be true for the conductance, because the edge brings in dependence 
on the integer part of a xy . A first example of this is the simple fact that the mean of the 
Hall conductance, that can be defined in a multiprobe geometry (with edges), depends on 
which transition is being studied, thus violating universality to this extent. The same may 
be true of the critical conductance fluctuations in geometries with edges. On the other 
hand, for a Corbino disk, which has no edges, there should be full universality among the 
integer quantum Hall transitions. Clearly, it would be of interest to study this numerically 
or experimentally. One way to do so numerically would be using the Chalker-Coddington 
model P3H , with additional co-moving edge channels coupling to the edges by hopping terms 



to obtain |cr° | > 1. 

Returning to the perturbative, metallic regime, a xx ^> 1, we expect that similar phenom- 
ena to those studied here in 2D should occur also in higher dimensions, for example in 3D. 
No isotropic topological term containing only two gradients is possible in higher dimensions. 
However, the Hall conductivity should make an appearance in the NLaM effective action, 
since it is a part of the measurable conductivity. It appears in a generalization of the 2D 



action to 3D 58 



S = - l -a° J d^rtild^Qd.Q} - ±a° B J d 3 re x , u n x tr [Qd,Qd u Q\. (6.1) 

Here n is a unit vector in the direction of the magnetic field B, and cr° and a° H are the 
diagonal (dissipative) and Hall conductivities, respectively. Thus the action is anisotropic, 
because the B field specifies a direction. (However, for simplicity we neglected the possible 
anisotropy in the diagonal conductivity a .) The action can be viewed as resulting directly 
from considering layers stacked perpendicular to the magnetic field, each of which has a 
Hall conductivity and is descibed by the 2D NLaM action, plus a transition amplitude for 
electrons hopping between the layers. Such models have recently been studied numerically 
||59|| . For systems with boundaries, the <j° h term in the action now leads in perturbation 
theory to phenomena similar to those in 2D, such as a tilted boundary condition, a depen- 
dence of (g) and the conductance fluctuations on a^/a , and so on. Thus in 3D, and also in 
still higher dimensions, the conductance fluctuations in general depend on the Hall ratio (or 
angle). However, for the localization transition in 3D, which would be expected to be in the 
unitary class since time-reversal symmetry is broken by the magnetic field, we suspect that 
the a° H term is irrelevant at the critical fixed point, so that the properties of the transition 
are universal, independent of the bare Hall ratio, at least to the same extent as in 2D, as 
discussed above. Similarly to 2D, the a° H term contributes to the action of configurations in 
which each layer has a non-zero instanton number (insofar as this number is well-defined, 
if the system has boundaries). In 3D, there also exist topologically-stable point-singular 
configurations of the Q field (known as "hedgehogs" in the literature), which may be viewed 
as points at which the instanton number changes from one layer to the next. The <j° h term 
counts the number of layers with each value of the instanton number, and thus is sensitive 
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to the presence and location of the hedgehogs. However, in the case of NLaM's studied in 
connection with antiferromagnets, it appears that the hedgehogs are irrelevant as far as the 
critical properties are concerned, even though they may affect the behavior in the phases on 
either side of the transition (see, e.g., Ref. [pO] ). Therefore, we suspect that, while the a° H 
term plays a role in the metallic phase, and also (after renormalization) in the 3D quantized 
Hall phase of layered systems f59j, it may have no effect on the critical properties, except 
perhaps for the conductance when edges are present. Clearly, these are questions that may 
repay further study. 
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APPENDIX A: BOUNDARY CONDITION AT THE HARD WALLS FOR THE 

WHITE-NOISE MODEL 



In this appendix we briefly derive the boundary condition on the diffusion propagator 
within the SCBA, by using current conservation at the reflecting walls. 

Using the SCBA equation (E — Hq — X^r))^? (r, r') = 5(r — r'), we can show that 



V- J + ~(r,r') = V- 



h 



G (r',r) 



-ifi 



2m P 



D G (r, r') 



-5(r-r')/u + G + (r,r')G (r',r) 



(Al) 



Let us define the G + G ladder diagram with n impurity lines as S + ' and the ladder 
diagram with one current vertex attached to the left as 

v L ' (n) (r,r') = u J rf 2 rx J + "(r,r 1 )5 + '-' {n) (r 1 ,r'). 



Using the above property of J + , we get the following recursive relation: 

V ■ v i,(n) = r [<?+ - S + ( n+1 )" . 



Summing up all ladder diagrams, we get 

V- v L (r,r') 



-5{r - r') 



(A2) 



(A3) 



where v L (r,r') represents the J + ~S + '~(r, r'). This shows that, on the finest length-scale 
resolution, (<r^(r, r'))scBA obeys V • <r(r, r') = 0, for r sufficently far from r'. 

For r at the reflecting boundary, the normal component of j outside the sample is zero. In 
the presence of a boundary current, which, from a coarse-grained, large-scale point of view, 
can be treated as (^-functions 5(y — W), 5{y) in the components tangential to the edge, the 
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surface integral of the current emerging from a small box centered on the top edge reduces 
to 



/•W+0+ 

/ dyd x (a x „(v, r')) SC BA - (<v( r > r '))scBA| y =w = 0, (A4) 

JW-0+ 

for r 7^ r'. Thus any normal current (just inside the edge) must be converted to a 5-function 
tangential current at the edge. This condition was discussed for the current j' in Sec. IA. 
Within the SCBA, it leads (using ( |2.20|) ) to the conclusion that it is 7 = cr® /a® x , not 
a ly l a xxi which appears in the boundary condition (P-.16Q on the diffusion propagator d. A 



similar argument holds for the r' dependence. The extension of this discussion to include 
the situation r = r' is given in Sec. IIIC. 

APPENDIX B: PROOF OF THE CURRENT CONSERVATION IDENTITIES 

WITHIN SCBA 

We will show that within SCBA, <7 MJ ,(r, r') satisfies the constraints imposed by current 
conservation. We start with the a + ~ term. We can write the ladder diagrams in the following 
fashion, 

h A r 

^ladder ( r , r ') = ^« J d 2 r 2 v L (r,r 2 )J t+ (r',r 2 ). (Bl) 

Using the recursive relation ( |A2|) for v L '^ n \ and denoting S ,+_ '^ n - ) as v R '( n \ we get another 
recursive relation: 

V ■ a+-> (n) (r, r') = ^(£+ - E^v^^r, r') - v*' (n) (r, r')]. 

Summing up all the ladder diagrams, we get 



V ■ ^ + -(r, r')>EE fefo r % P N '(r>, r) - P N (r, r') D G^r', r) 

n N N' 

where P^(r,r') is the projection operator onto the Nth. Landau level. The right hand side 
is a short-ranged function of |r — r'| , which we can treat as a 5-function. We can write 

V-cr + "(r,r / ) = c +- 5(r - r'), (B2) 

where 

n n < 

where v nn i is the matrix elements of the velocity operator, and 



for B = 0; 







^ for B ^ 0. 

To evaluate a ++ (r, r') and a (r, r'), we use the following trick: 
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„"(r, r<) = £ iET /(#) Ei Jim e , ^(r, r<; Bl , E 2 ), 
a-(r, r<) = £ ^ /(£<) ta^ ^"(r, ✓ I ft, ft), 



(B3) 



where a aa {E 1 ,E 2 ) involves ladder sum S aa (a = +, — ). Define V - R,aa '( n ) = S aa ^J aa , we can 
show that for the nth ladder diagram, 



2m, 



^E+^-E+^Jv 



fl,++,(n-l), 



ft* 

2m e 
"IF 

+5 £ ^i-^]v^ ++ 'W(r,r'). 

ft. 



r, r 



[E + (^ 1 )-E+(£; 2 )]v ii ' ++ 'W(r,r / ) 



(B4) 



One can see that there is cancellation between the nth ladder diagram and the (n + l)th 
ladder diagram. Similar relations can be derived for V • a >( n )(r, r'; Ei, E 2 ). Summing up 
all the ladder diagrams, taking the derivative over energy and then the limit E\ = E 2 — > E' , 
we can show that 

V ■ [cr ++ (r, r') + o (r, r')] = — c + ~5(r — r') 

" Sri dE 'f( E ')\v ++ > R (ry,E') -y-' R (r,v',E')]. (B5) 
We can see that V-<x + ~(r, r') is canceled by contributions from V-o~ ++ (r, r') and V-er (r, r'), 

h 2 f E 



2m, 



f dE'f(E')[v ++ ' R (r,v',E')-v~' R (v,r',E')} 



(B6) 



Since 



and 



we get 



V'-v fl - ++ 'W(r,r') 



2m. 



h 

2m 



r, r 



1 -[S ++ ^ n - 1 \r,i f ) - 5 ++ '^(r,r')] = 0, 
-[d~ > in - l) (r,r') -rf— • (n - 1) (r,r / )] = 0, 



V( ( r(r,r / ))scBA- V'= 0. 



(B7) 



Using eqs. B6 and B7 and the asymptotic property of the Green's function [42 



we can show finally that 



dS'L^ooGH^r') D G ± (r',r) = 0, 



V- |(a(r,r / ))scBA-rfS / = 0. 



(B8) 
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APPENDIX C: REMARKS ON EDGE STATES AND QUANTIZATION 



Here we return to the topological considerations of Sec. IIIB, and relate them to edge 
states and the quantization of the Hall conductance in the localized regime. The topological 
considerations of Sec. IIIB are closely related to the problem of setting up a path integral 
for a quantum spin, by which we mean an irreducible representation of the symmetry group, 
which is SU(2n) here (for a review, see e.g. |60U), and this connection is also utilised in the 
mapping from the Chalker-Coddington model (representing a network of edge states) to a 
quantum spin chain or the NLcrM |49[] (the connection between the latter two problems, 
and the relation to the quantum Hall effect, was discussed earlier [^TJ). In the quantum 
spin problem, we would take imaginary time, with a periodic boundary condition in the 
time direction, and the action would contain only the cr^'° terms from 5 [A], eq. (|3.8|); the 
system would be taken to be a disk, with the single edge corresponding to the world line of 
the quantum spin with its periodic boundary condition. For the two-probe geometry, this 
corresponds to regarding x as imaginary time, and the two edge channels are then a pair 
of quantum spins, with the spins fixed at Q = A at the initial and final "times" x — 0, L. 
In the absence of the rest of the action, quantum-mechanical consistency requires in either 
geometry that the coefficient a^y be quantized to integer values, for reasons closely related 
to the properties of "large" (topologically non-trivial) gauge transformations; for the case of 
SU(2), this corresponds to 2S = integer, as usual. Essentially, the argument says that, since 
the only degree of freedom in the problem is the value of Q on the edge, then its continuation 
into the interior, needed to write the topological term, is arbitrary, and the path integral 
should be invariant under a change in Q in the interior that does not affect the edge; such 
changes are the "large" gauge transformations. Since the change in the action under such a 
change is 2iria I J y , ° q for some integer q, this implies that cr^'° is an integer. This is related 



to arguments for quantization of the Hall conductivity, once localization sets in |29|,|3(J. In 
this case, we may imagine that the localized system is described by the NLcrM but with 
cr° x replaced by a renormalized value a xx equal to zero because of localization. Then a 
similar argument requires that the renormalized a]L is quantized to integer values. Thus 



quantization of the Hall conductance and quantization of spin are closely connected ^2 
This argument is also connected |30| with the gauge-invariance argument for quantization 
TOfl . The edge states, that are the only degrees of freedom able to transport current over 



large distances in the localized regime, correspond to the quantum spin (in the n — > limit). 
We note that from this point of view of the edge states, in which x plays the role of imaginary 
time, A x plays the role of an external magnetic field, in the sense of the familiar Zeeman 
coupling in the SU(2) case. It is coupled to Q, which corresponds to the spin operator, or 
the current operator for the edge state. Q corresponds for SU(2n) to the three-component 
unit vector f2 that describes an SU(2) spin, which can be obtained explicitly by writing, for 
n = 1, Q = Q ■ r, where r is the vector of Pauli matrices. For SU(2), the corresponding 
operators in the quantum theory, after rescaling to absorb the coefficient analogous to our 
cr^y , are the familiar operators S, which generate SU(2) rotations and are conserved when 
the Hamiltonian is SU(2) invariant. In the presence of the vector potential A, which enters 
multiplied by the magnitude of the spin, 5*, to give the Zeeman coupling, just as in our action, 
the equation of motion (or Ward identity) for a single quantum spin describes the familiar 
precessional dynamics, which can therefore also be viewed as the covariant conservation of 
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the spin. 

By contrast, for the full action in the weak coupling regime a xx ^> 1, where the rest of 
the action depends on the form of Q in the interior, there is no reason why either a xy or 
a^y should be quantized, in accordance with our physical expectations. The same applies 
to derivations of S[A] starting from the network model ||49|| , which also represents only the 
Fermi energy response, except that in this case the splitting of cr° into cr^° and a^l' is 
a matter of an arbitrary definition, as we discussed for the linearized model in Sec. IB. In 
the network model, the links of the lattice can be viewed as quantized edge channels, and 
these are the only degrees of freedom, so the coupling of A to these links is of the edge form 
discussed above, but with a*y'° replaced by the quantized value 1. Of course, the coarse- 
grained values of a xx , a x 'y, and <t^'° are determined by the parameters of the vertices in 
the network model, and by the definition of the coarse-grained currents [0 , so they are not 
quantized. 



APPENDIX D: COMPUTATION OF THE VARIANCE DIAGRAMS 

The conventional diagrams Figs. |]a and b depend on 7 through the diffusion propagator. 
Using the boundary perturbation expansion, we get 

Tr(dd) = Tr{d°d°) + Tr{d° d°Bd°) + Tr{d°Bd° d°) 

+ 2Tr(d° d°Bd°Bd°) + Tr{d°Bd° d°Bd°) + 0( 7 3 ), (Dl) 
Tr(dd T ) = Ti(d°d°) + Tr(d° d°B T d°) + Tr(d°B T d° d°) 

+ 2Tr(rf° d°B T d°B T d°) + Tr(d°B T d° d°B T d°) + 0( 7 3 ), (D2) 

where the matrix B has the following elements in the basis of 0° m : 

v _ 7 8mm' 

{Tim LD 77,777.) — ——— — ~0 m + m i o dd0n+ri, even- 
lY L {m') z — m z 

The linear term in 7 is zero, because the matrix B is anti-symmetric. We get 



1 64 

(^ = -[6/ 2 -7 2 -( 12 /3-2/ 4 )], (D3) 



where 



h iw) 55 (m 2 + n 2 L 2 /lF 2 ) 2 ' (D4) 
f f—\ ST ST 1 1 (mm 1 ) 2 

h J^J^ (m 2 + n 2 ^f m> 2 + n> 2 " \{m>) 2 - m ^™'^^— 

(D5) 

( L \ sr^ sr^ * * (mm') 2 

2] 2 Om+m' .oddOn+n' ,e 

(D6) 



•Z" 4 ( Ti7 ) / o , n r2 in , ,n , ,n r2 1, \i ,\o 91 o "m+m' ,odd^n+n' ,even 

V^/ n4^=in'=5X'=i (m 2 + n 2 -^) 2 (m' 2 + n' 2 ^) 2 [(m') 2 -m 2 ] 2 
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The upper bounds for n and m in all sums are n max ~ W/l and Tn, max ~ L /I. For a square 
sample with L = W, we have f 2 — 1.51, 

(59)l b ,L=w = At 9 - 06 - 2 - 35 ^] + 0(7 3 )- (D7) 
Diagrams ^e and f give 

^lf = 2 ^ 2 W S E E ^ 

ni ,mi=oaa 712, m 2 =odd n 3 =U,m 3 =l x ^ 
1 1 1 

X m 2 + nfL 2 /iy 2 ml + 77 2 L 2 /PF 2 (ml + 77 2 L 2 /W 2 ) 2 
x {3 (m 2 3 D 3 D 2 + n 2 3 L 2 /W 2 £>i£> 4 ) + 3 (m x m 2 D^D 2 + ni n 2 L 2 /W 2 D 1 D ( 

+ 10 (m im3 D 7 D 2 + nmsLV^DxlV) } , (D8) 

where 

Di(mi,m 2 ,m 3 ,m 3 ) = 5 mi ,m 2 — — ^mi,2m 3 +m 2 ~~ -fim 1 ,-2m a +m 2 + " ^"1 >2m3 -m2 > 
D 2 (rii,n 2 , 713,71%) = 5 nit n 2 + -^ni,2n 3 +n 2 + ^n 1 -2n i +n 2 + ^ ^l,2n 3 -n 2 ) 

D3(lTli, m 2 , 777,3, 777.3) = ^mi,m 2 + 2^ m i' 2m 3+ m 2 _ 2^ mi > _2m3+m2 _ T^" 11 ' 2 " 13- " 12 ' 

111 

D^fli, U 2 , U 3 , 77,3) = 5 nii „ 2 — -<5rM,2n 3 +n 2 + ~ ,-2« 3 +n 2 — ~ ^ ra i >2n 3 -n 2 ) 

1 1 1 

D^ifTll, 7772, m.3, m.3) $mi,m 2 2 ^Wi ,2m 3 +m 2 ^ ^'^li - 2m 3 +m 2 2 i2m 3 — m 2 > 

111 

jD 6 (n!, 77-2, 77-3, 77 3 ) = 5„ ljn2 + -^ni,2n 3 +n 2 + 7^n 1 ,-2nz+n 2 ~ ^n 1 ,2n i -n 2 -, 

111 

Df{m\, 777.2, W3, ?773) — ~ 2^ m l' 2m 3+ m 2 + 2^ m l' _2m 3+ m 2 + 2^ m l' 2m 3 _m 2! 

Ill 

Ds(ni, n 2 , 773, 773) = -<5 nii 2n 3 +n 2 ~ -^ n i ~ 2ri 3+ n 2 + 2^ ni > 2n3 ~ U2 ' (D9) 

This term has a logarithmic part (it diverges with system size as log(L/Z)). It comes from 
the first term in the curly brackets, when the two derivatives of the 4-point interaction 
are applied to the closed loop of two diffusion propagators. The second term results from 
applying the two derivatives to the two external propagators. The third term arises when 
one of the derivative is applied to the closed loop, one is applied to the external propagator. 

Diagram Figs. |]g and h are both logarithmic. They are of opposite signs, but the 
amplitude of diagram f|g, which is positive, is twice that of |h. We get 

^9)l h = ^i 2 ^- 2 E E E ^Jr 

" vr ni,mi=oddn2,m2=oddn 3 =0,m3=l 1 2 

111 , 
* 2 I Z2TT7^—2 I „n a / m „2 I .2r,/„^ l^. (D10) 



m 



\ + n\L 2 jW 2 m 2 + n 2 L 2 /W 2 m 2 + n\L 2 /W 2 
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Diagram |]j and IJj' also have logarithmic divergence. We get from diagram f|j 

(<H 2 = fr7 2 ^ E E E E ^r^r 

ni,mi=oddn2,m2=oddn:j=0,m i =ln4=0,m4=l 1 A 

1 1 m 3 m 4 

X m\ + nlL 2 /W 2 m\ + n 2 2 L 2 /W 2 ml + n 2 3 L 2 /W 2 m\ + n\L 2 /W 2 

^ [^mi ,m4±ni3 ^mi,m3- 7714] [^m2,m3±m4 ^m2,ni4- 7713] 

(Dll) 

We get from Uf 

( ^ 2, = "^ 72 ^ E E E E 

ni,mi=odd n,2, m,2=odd 113=0, m,3=l n4=0,m4=l 1 z 

1 1 m\ 1 

X mf + nfL 2 / W 2 m| + n 2 2 L 2 /W 2 ml + n§L 2 /W 2 m\ + n^L 2 /W 2 

^[^mi,iri4±m3 ^mi ,013— ^14] [<^m2 ,m4±m3 $1712, 7713—7714] 

^ [^ni,n4±n3 "I - ^ni ,773—714] [^77,2 ,rt4±7i3 "I - ^77,2,713 — 714] • (1)1 2 J 

Both diagrams U and j' are negative. Their logarithmic parts combine to cancel those from 
diagram e, f, g and h. The variance, which is the sum of a-j', is finite. 
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FIGURES 



FIG. 1. The two-probe geometry. 

FIG. 2. The Chalker-Coddington network model. Each unit cell contains four distinct links, 
A, B, C, D. The tilted boundary condition arises from the fact that along each link, the random 
walk is along only one direction. 

FIG. 3. (a) The SCBA single-particle Green's function. It sums up all the non-crossing dia- 
grams. The thin line denotes G°, the Green's function in the absence of disorder. The thick line 
denotes the SCBA Green's function G. (b) The ladder sum for the SCBA two-particle Green's 
function, (c) The diagrams for the SCBA bilocal conductivity tensor. 

FIG. 4. The diagrams for the variance of conductance to leading order in 1/(T® X and to order j 2 . 
The shaded polygons are vertices. The lines connecting the vertices are the diffusion propagators. 

FIG. 5. The equivalence between the diagrammatic approach and the NLctM approach. One 
vertex from the NLctM, 5 2 S/6A£~[ 1 5A~ jj, is equal to the sum of four diagrams in the diagrammatic 
approach. 
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